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PREFACE 
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B.  Barclay,  and  J.  W.  Workman  of  the  Agbabian  Associates  staff.  Technical 
monitor  of  the  contract  for  the  Department  of  Defense  Explosives  Safety  Board 
is  Dr.  T.  A.  Zaker. 
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SUMMARY 

This  rerort  describes  Che  results  of  dynamic  response  analyses  of 
magazine  headwalls  in  three  tests,  namely,  ESKIMO  I,  It,  and  IV.  The  finite 
element  models  of  the  headwalls  were  subjected  to  blast  loadings  selected  on 
the  basis  of  measured  data  and  theoretical  investigations.  The  responses  of 
the  headwalls  were  obtained  using  the  INSLAB  code,  a nonlinear,  dynamic 
finite  element  computer  program.  The  behavior  of  the  headwall  material  was 
described  by  a bilinear  moment-curvature  relationship,  and  the  supporting  soil 
was  simulated  by  a series  of  nonlinear  springs  and  linear  dampers.  The  results 
of  the  calculations  were  compared  with  available  test  data.  In  general,  the 
predicted  behavior  of  the  headwalls  agreed  very  well  with  their  observed 
behavior.  The  present  analysis  of  the  south  headwall  in  ESKIMO  I,  using  the 
refined  finite  element  showed  better  agreement  with  the  test  results 

than  those  from  the  previous  analysis,  which  used  a coarse  mesh.  The  results 
for  the  east  headwall,  ESKIMO  II,  showed  that  the  computed  permanent  displace- 
ment contours  are  quite  similar  in  shape  to  Chose  measured.  The  results  for 
the  northeast  headwall,  ESKIMO  II,  indicated  that  the  concrete  beams  around 
the  doorway  were  effective  in  reducing  the  maximum  displacement  of  the  headwall 
and  that  the  biparting  and  sliding  type  of  door  was  superior  to  the  double- leaf 
and  hinged  type  for  resisting  blast  loads.  The  ccmputea  motions  of  the  north- 
east headwall,  ESKIMO  IV,  compared  very  well  with  the  corresponding  test  data. 
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SECTION  1 

INTRODUCTION  AND  BACKGROUND 

This  report  covers  a part  of  a general  program  sponsored  by  the 
Department  of  Defense  Explosives  Safety  Board  (DDESB)  to  determine  safe  inter- 
magazine separation  distances  for  various  orientations  of  magazines  for  storing 
chemical  explosives.  The  main  objective  of  the  program  is  to  determine  the 
minimum  intermagazine  separation  distances  so  that  an  explosion  in  one  magazine 
(donor)  will  not  cause  an  explosion  in  an  adjacent  magazine  (acceptor).  This 
report  describes  a study  to  predict  analytically  the  behavior  of  magazine 
headwalls  and  to  compare  analytical  and  test  data. 


Preliminary  tests  had  indicated  that  specifications  for  minimum 
separation  distances  between  magazines  were  excessively  conservative  for  some 
orientations.  Increasing  land  costs  and  siting  problems  made  it  desirable  to 
justify  reducing  the  rear-to-f ront  separation  distance.  Ea.  ier  tests  had 
indicated  that  earth-covered,  steel-arch  igloo  magazines  can  be  safely  spaced 
side  by  side  at  a distance  in  feet  of  1.25  W1y^,  in  which  W is  the  weight  in 
pounds  of  the  high  explosives  in  storage.  On  the  other  hand,  very  few  test 
data  had  been  developed  foF  determining  the  minimum  safe  distances  for  other 
magazine  orientations. 

Starting  in  1971,  the  DDESB  sponsored  a series  of  large-scale  maga- 
zine explosion  experiments,  called  the  ESKIMO  series,  for  the  purpose  of  estab- 
lishing  minimum  separation  distances  for  earth-covered,  steel-arch  magazines. 

To  date,  four  tests  designated  as  ESKIMO  I,  ESKIMO  II,  ESKIMO  Ml,  AND 
ESKIMO  IV  have  been  conducted.  These  tests  are  described  briefly  in  the 
following  sections. 


1.1  DESCRIPTION  OF  ESKIMO  I TEST 

ESKIMO  I was  a large-scale  magazine  explosion  experiment  conducted 
on  8 December  1971  at  the  Randsburg  Wash  Test  Range  of  the  Naval  Weapons  Center, 
China  Lake,  California.  Four  earth-covered,  steel-arch  magazines  were  exposed 
to  the  explosion  of  the  contents  of  a similar  magazine.  The  donor  magazine 
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contained  200,000  lb  of  high  explosives.  The  four  acceptor  igloos  faced  the 
donor  and  were  located  at  various  distances  ranging  from  73  ft  to  1 61  ft,  as 
shown  in  F’gure  1-1.  The  distances  73  ft  and  1 6 1 ft  correspond  to  1.25  W1^ 
and  2.75  W1^,  respectively,  where  W is  the  weight  in  pounds  of  the  explo* 
sive  in  the  donor  magazine.  Two  concrete  block  structures  simulating  cne 
particular  type  of  Air  Force  aboveground  magazine  were  also  placed  in  the 
area  at  distances  of  117  ft  (2  from  the  donor  igloo,  as  shown  in 

Figure  1-1. 

The  four  acceptor  steel-arch  igloos  were  built  in  accordance  with 
the  standard  drawing  shown  in  Figure  1-2.  Each  of  the  igloos  was  25  ft  wide 
by  14  ft  high,  with  the  length  limited  to  20  ft.  Steel  wing  walls  were  used 
in  the  test  in  lieu  of  concrete.  The  igloos  were  covered  by  a 90%-compacted 
earth  surcharge  to  a depth  of  2 ft  at  the  top  of  the  arch. 

The  primary  objective  of  this  test  was  to  establish  a safe,  practi- 
cable minimum-separation  distance  between  steel-arch  magazines  for  face-on 
exposures.  Another  major  objective  of  the  test  was  to  determine  the  fragment 
and  debris  hazard  from  mass  detonatiop  of  explosive-filled  projectiles  in  an 
earth-covered  magazine. 

• 

The  damage  to  structures  varied  from  minor  headwal 1 damage  at  the 
south  igloo  to  complete  destruction  of  the  concrete  block  structures.  The 
status  of  the  acceptor  charges  after  the  test  indicated  a range  from  no 
explosion  or  burning  in  the  south  igloo  to  complete  explosion  or  detonation 
of  all  charges  in  the  east  igloo. 

Permanent  deflections  of  the  order  of  several  inches,  accompanied 
by  yield-line  formation,  have  been  noted  in  some  of  the  surviving  headwal Is. 

On  the  other  hand,  photographic  evidence  indicates  that  the  steel  plate  doors 
in  r o igloos  were  driven  in  with  considerable  velocity  before  coming  to  rest, 
while  remaining  partially  attached  to  frames.  In  the  forward  exposure,  com- 
plete failure  of  the  hinges  occurred,  and  the  doors  were  driven  violently 
against  the  rear  wall  of  that  magazine  (20-ft  depth,  wall  to  wall). 
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Refsrence  1 includes  e detailed  description  end  the  conclusions  of 
the  E'.KIMC  I test.  Test  results  indicated  that  formerly  applicable  separation 
distances  could  be  significantly  reduced. 

Agbabian  Associates  (AA)  performed  dynamic  response  analyses  of  the 
headwatls  in  the  ESKIMO  I test,  using  a finite  element  model  of  the  hecdwall? 
(Ref.  2).  The  analyses  were  performed  with  the  INSLAB  code,  a nonlinear  finite 
element  computer  program  (Ref.  3)*  The  objective  of  the  analyses  was  to 
understand  the  behavior  of  the  reinforced  concrete  headwalts  subjected  to 
blast  loading  and  to  use  this  knowledge  in  reducing  reliance  on  full-scale 
proof  testing. 

1.2  DESCRIPTION  OF  ESKIMO  II  TEST 

■The  ESKIMO  II  test  was  conducted  on  22  May  1973  at  the  Randsburg 
Wash  Test  Range  of  the  Naval  Weapons  Center,  China  Lake,  California.  Five 
earth-covered,  steel-arch  magazines  were  exposed  face-on  all  at  the  same 
distance  from  the  explosion  source,  as  shown  in  Figure  1-3.  The  igloos  to  the 
north,  south,  and  west  were  those  remaining  from  the  ESKIMO  I test.  However, 
the  igloos  to  the  south  and  west  were  fitted  with  new  door  designs.  All  the 
magazine  structures  were  earth-covered,  semicircular  corrugated  steel  arches 
except  the  northeast  igloo,  which  was  constructed  with  a new  type  of  noncir- 
cular steel  arch.  Table  1-1  shows  details  of  all  igloos  in  the  ESKIMO  II  test. 

The  explosion  source  consisted  of  72  tri tonal -fi 1 led,  750-lb  bombs 
in  two  triangular  stacks,  in  base-to-base  contact.  This  source,  containing 
27,800  lb  of  explosives,  had  an  expected  TNT  equivalent  weight  of  approximately 

24.000  lb.  The  source  was  designed  to  produce  an  impulse  load  of  1100  psi  msec 
on  the  headwalls  of  the  test  magazines  located  at  147  ft  from  the  source.  This 
was  the  value  of  impulse  estimated  to  act  on  the  lowest  third  of  an  acceptor 
igloo  from  explosion  of  an  earth-covered  magazine  filled  to  capacity  (i.e., 

500.000  lb  net  explosive  weight  permitted  in  one  igloo  by  the  current  standards) 
at  a rear-to-f ront  separation  distance  of  2.0  x-W^^,  or  159  ft,  where  W is 
the  weight  (lb)  of  charge  in  the  magazine  filled  to  capacity.  Subsequently, 
however,  free-field  measurements  made  to  the  rear  of  the  ESKIMO  III  donor 
magazine,  and  loadings  observed  in  ESKIMO  IV,  have  indicated  an  impulse  load 
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of  about  1200  psi  msec  from  only  350,000  1b  detonated  in  an  igloo  at  this 
scaled  distance.  In  ESKIMO  IV,  this  level  of  loading  was  produced  1 4 7 ft 
from  a bare  TNT  hemisphere  weighing  37,000  1b. 

• 

Piezoelectric  and  strain-gage  pressure  transducers  were  positioned 
in  and  near  the  headwalls  of  the  Zest  magazines.  Accelerometers  and  linear 
motion  transducers  were  mounted  on  the  doors  and  headwalls  of  the  structures. 
Self-recording  mechanical  pressure  gages  were  placed  at  distant  stations. 

Door  motions  were  observed  by  a high-speed  camera  in  the  interior  of  each 
of  the  magazines. 

Because  the  near-field  blast  loading  exceeded  that  planned,  the 
igloo  structures  were  subjected  to  an  overtest.  Despite  this  overtest,  the 
large,  single-leaf  sliding  door  on  Igloo  D (south)  withstood  the  blast  loading 
without  breakup  or  severe  distortions.  Likewise,  the  Stradl ey-type  headwall 
used  for  Igloo  B (northeast)  incurred  only  an  acceptable  degree  of  damage. 

The  test  reaffirmed  the  need  for  achieving  a balance  in  the  strength  of 
headwalls  and  doors. 

Reference  4 includes  a detailed  description  of  the  ESKIMO  II  test. 

1.3  DESCRIPTION  OF  ESKIMO  III  TEST 

The  ESKIMO  III  test  was  conducted  on  12  Jun*-l-9.74  at  the  Randsburg 
Wash  Test  Range  of  the  Naval  Weapons  Center,  China  Lake,  California.  Five 
earth-covered,  steel-arch  magazines  were  exposed  in  the  test  to  the  explosion 
of  the  contents  of  a similar  magazine,  as  shown  in  Figure  1-4.  Igloos  B, 

D,  and  E were  those  remaining  from  the  ESKIMO  II  test.  Igloo  A and  the  donor 
magazine  were  new  structures.  Igloo  A structure  substituted  a less  expensive, 
light-gage,  deep-corrugated  arch  in  place  of  heavy  structural  plate;  other- 
wise, it  was  identical  to  the  standard  steel-arch  igloo.  This  igloo  would 
furnish  a test  of  the  opinion  that  the  earth  cover  is  the  most  important 
factor  in  preventing  explosion  communication.  In  addition,  this  igloo  was 
intended  to  provide  a direct  qualification  test  of  the  lightweight  igloo. 

The  donor  magazine  was  built  to  the  same  specifications  as  Igloo  A.  Table  1-2 
shows  details  of  all  igloos  in  the  ESKIMO  III  test. 
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The  explosion  source  consisted  of  350,000  lb  of  tri tonal  in 
750- 1 b bombs  in  the  donor  igloo  shown  in  Figure  1-4. 

Incident  peak  pressure  was  measured  at  the  ground  surface  in  front 
of  each  acceptor  igloo.  Reflected  peak  pressure  was  measured  from  face-on 
blast  gages  set  in  each  headwall.  A row  of  surface  gages  measured  the 
incident  peak  pressure  over  the  earth  on  the  oval  arch  igloo  and  on  the  light- 
gage  circular  arch  magazine.  Acceleration  and  displacement  measurements  were 
made  on  the  arches  of  these  igloos  so  that  should  they  fail,  the  rapidity  of 
failure  could  be  used  to  estimate  whether  there  would  have  been  a risk  of 
explosion  communication  to  the  contents.  Detailed  damage  surveys  and  permanent 
displacement  measurements  were  made  after  the  test. 

The  main  objective  of  the  ESKIMO  III  test  was  to  expose  an  earth- 
covered,  noncircular,  corrugated  steel  arch  magazine  (Igloo  B)  and  a new 
circular,  light-gage,  deep-corrugated  steel  arch  magazine  (Igloo  A)  to  explo- 
sion of  an  adjacent  magazine  at  the  minimum  side-to-side  spacing  now  permitted 
by  standards.  In  addition,  magazine  headwall  and  door  structures  were  tested 
at  several  other  distances  arid  orientations  of  interest. 

A preliminary  inspection  of  the  structures  after  the  test  snowed 
that  the  doors  of  Igloos  A,  B,  D,  and  E were  all  blown  open.  The  headwall  of 
Igloo  E experienced  the  most  severe  damage,  and  permanent  inward  displacements 
were  as  much  as  1.5  ft.  Except  for  local  damage  near  one  edge  of  the  door 
opening,  the  headwall  of  Igloo  D experienced  generally  slight  damage.  The 
remaining  headwalls  showed  only  minor  damage.  Additional  details  on  the  test 
are  provided  in  Reference  5. 

1 .4  DESCRIPTION  OF  ESKIMO  IV  TEST 

The  ESKIMO  IV  test  was  conducted  on  10  September  1975,  also  at  the 
Randsburg  Wash  Test  Range,  Naval  Weapons  Center,  China  Lake,  CaUfornia.  Its 
principal  objective  was  to  demonstrate  the  resistance  of  a newly  designed 
headwa 1 1 -and-door  combination  to  explosive  blast  loading  from  detonation  of 
the  contents  of  another  magazine. 
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The  ESKIMO  IV  test  array  consisted  of  three  magazine  structures 
each  facing  the  explosion  source  147  ft  away,  as  shown  in  Figure  1-5*  The 
primary  target  structure  was  Igloo  B to  the  northeast,  consisting  of  a single- 
leaf sliding  door,  horizontally  spanned,  mounted  on  a modified  headwall  of  the 
standard  Stradlfey  magazine. 

Igloo  0 (east)  was  rebuilt  with  the  two-leaf,  hinged,  steel-plate 
door  and  the  headwall  of  the  standard  circular  steel  arch  magazine.  Exposed 
to  the  same  level  of  loading  as  the  primary  target,  it  served  as  a control 
structure  to  demonstrate  directly  the  relative  strength  of  the  primary  target 
(Igloo  B) . 


Igloo  E to  the  west  was  rebuilt  with  the  headwall  of  the  standard 
circular  steel  arch  magazine,  but  fitted  with  a single-leaf  sliding  door. 

This  combination  was  v»«ted  inconclusively  in  ESKIMO  II,  but  the  test  had 
indicated  a serious  Imbalance  in  strength  between  the  door  and  headwall. 

The  test  utilized  a nearly  ideal  explosion  source  to  generate  blast 
* loadings.  It  afforded  the  opportunity  for  more  extensive  source  diagnostics 
and  dynamic  response  measurements  on  the  target  structures  than  in  previous 
tests  of  the  ESKIMO  series.  This  in  turn  would  enable  an  extensive  correla- 
tion  between  measured  and  calculated  responses. 

The  explosion  source  consisted  of  8— 1 b rectangular  blocks  of  TNT 
stacked  in  the  form  of  a hemisphere,  containing  a total  of  37,000  lb  of  TNT. 
The  explosives  were  left  exposed,  rather  than  surrounded  by  an  earth  cover 
or  barriers. 

Pressure  gages  were  installed  in  the  headwalls  and  in  the  ground 
nearby  to  measure  blast  pressures  on  the  »headwa  1 1 s . Accelerometers  and  linear 
displacement  transducers  (LVDTs)  were  used  to  measure  structure-response  time 
histories  of  the  headwalls.  Sel f-recording  pressure  gages  were  installed  to 
meacure  pressure  loading  at  the  far-field  stations.  Deta'ls  on  the  locations 
of  the  pressure  gages,  accelerometers,  and  LVDTs  are  discussed  in  Section  4. 
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NOTE:  IGLOOS  A,  C,  AND  F WERE 

OLD  STRUCTURES  FROM  THE  PREVIOUS 
ESKIMO  TESTS  THAT  WERE  HOT 
CONSIDERED  IN  ESKIMO  IV. 


FIGURE  1-5-  LAYOUT  OF  TEST  STRUCTURES  FOR  ESKIMO  IV  MAGAZINE  SEPARATION  TEST 
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Preliminary  results  (Ref.  6)  showed  that  the  doors  of  Igloos  0 
and  E were  blown  in.  The  headwall  of  Igloo  E sustained  extensive  damage. 

The  headwalls  of  Igloos  8 and  0 and  the  door  of  Igloo  B sustained  little  or 
no  damage. 

1 . 5 OBJECTIVES  AND  SCOPE 

The  overall  objective  of  the  present  study  was  to  develop  confi- 
dence in  analytical  techniques  for  predicting  the  dynamic  response  of  concrete 
headwalls  and  steel-door  structures  to  air-blast  overpressures  resulting 
from  explosion  of  the  contents  of  an  adjacent  magazine.  To  achieve  the  objec- 
tive, the  analytically  predicted  responses  of  the  headwalls  were  compared 
with  experimental  data  from  the  ESKIMO  tests.  The  following  tasks  were 
defined: 

Task  1 

Calculate  the  response  of  the  headwall  and  door  of  the  south  igloo 
in  the  ESKIMO  I test,  using  a finite  element  model  wiLn  improved  resolution. 

Task  2 

Cal  culate  the  response  of  the  headwalls  and  doors  of  the  east  (C)  and 
northeast  (B)  igloos  in  the  ESKIMO  II  test. 

Task  3 

Predict  response  of  the  strengthened  Stradleytype  headwall  and 
single-leaf  sliding  door  of  the  northeast  (B)  igloo  in  the  ESKIMO  IV  test.* 


-Because  of  tight  schedules,  this  task  could  not  be  completed  before  the  test. 
Nevertheless,  the  analysis  was  performed  before  any  results  from  the  test 

became  available. 
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SECTION  2 

PRESSURE  LOAOING  ON  HEADWALLS 

This  section  presents  air-blast  pressure  histories  used  in  the 
dynamic  response  calculations  of  headwall  and  door  systems.  Provisions 
for  measuring  air-blast  overpressure  data  on  headwalts  were  made  in  the  plans 
for  all  shots  in  the  ESKIMO  test  series.  However,  the  measured  data  were 
generally  meager  and  of  low  quality.  Therefore,  data  from  other  tests  and  from 
theoretical  considerations  were  used  to  design  the  pressure  histories  for 
the  analytical  calculations  of  the  ESKIMO  tests. 


2.1  PRESSURE  HISTORIES  FOR  ESKIMO  I TEST 

The  pressure  time  histories  for  the  south  igloo  in  ESKIMO  I used  in 
the  present  study  were  the  same  as  the  corresponding  histories  calculated  in 
the  previous  study  (Ref.  2).  These  histories  are  shown  in  Figure  2-1.  The 
three  pressure  histories  represent  loading  for  the  three  zones  of  the  headwall, 
as  shown  in  Figure  2-2,  and  differ  only  during  the  unloading  phase  of  the 
pressure  history.  The  differences  reflect  the  influence  of  unloading  signals 
propagating  downward  from  the  top  of  the  headwall.  The  use  of  just  three 
different  zones  was  considered  adequate  for  the  finite  element  calculations. 

All  the  pressure  histories  incorporate  a millisecond  (msec)  rise-time  ramp 
front  for  compatibility  with  the  integration  time  step. 

Additional  discussions  on  the  calculation  of  the  pressure  histories 
for  ESKIMO  I may  be  found  in  Reference  2. 
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FIGURE  2-1.  REFLECTED  OVERPRESSURE  ON  SOUTH  IGLOO,  ESKIHO  l 
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2.2  PRESSURE  HISTORIES  FOR  ESKIMO  II  TEST 

Although  redundant  air-blast  instrumentation  was  provided  at  all 
five  of  the  test  igloos  in  ESKIMO  II,  obtaining  reliable  detailed  descriptions 
of  the  blast  loads  applied  to  the  magazine  headwalls  was  difficult.  Records 
from  most  of  the  gages  were  reported  to  be  "difficult  to  read  and  considered 
unreliable."  Peak-incident  (side-on)  and  reflected  overpressures  measured 
at  the  east  and  northeast  igloos  were  consistent  with  other  experimental  data 
and  with  the  theory  of  gas  dynamics  for  shock  reflection.  On  the  basis  of 
the  measured  values  and  this  consistency,  a peak-loading  pulse  pressure  of 
260  psi  was  selected  for  the  ESKIMO  II  calculations. 

For  the  five  test  igloos  the  measured  impulse  varied  by  a factor  of 
two  from  the  lowest  to  the  highest  value.  The  nearly  twofold  variation  in 
measured  inpulse  may,  in  part,  oe  due  to  the  lack  of  rotational  symmetry  of 
the  charge  and  surrounding  revetment.  The  impulse  for  the  east  and  west 
igloos  is  well  below  that  for  the  other  igloos.  However,  significant  (20% 
to  30%)  north-south  and  east-west  variations  in  impulse  are  present.  The 
impulse  selected  for  the  headwall  loading  pulse  (1720  psi  msec)  is  considered 
representative  for  the  east  and  northeast  igloos. 

No  positive-phase  duration  data  were  reported  for  the  pressure 
histories  measured  at  the  headwalls.  Pressure  histories  were  not  available 
for  examination  to  establish  durations.  Scaling  of  previous  magazine-safety 
test  data,  general  high-explosives  and  nuclear-test  data,  and  computational 
data  suggested  durations  of  15  to  35  msec.  Consideration  of  typical  pressure 
history  pulse  shapes,  the  selected  peak  pressure,  and  the  desired  impulse  led 
to  the  selection  of  a 30-msec  duration.  The  resulting  pulse  developed  for 
loading  the  east  and  northeast  igloos  of  ESKIMO  II  is  shown  in  Figure  2-3. 

The  three-zone  description  of  the  pressure  loading  of  the  headwall 
was  retained  for  ESKIMO  II.  The  pressure  histo-ies  for  Zones  1 and  2,  which 
experience  earlier  unloading,  are  shown  in  Figure  2-3  superimposed  on  the 
pressure  history  for  Zone  3 derived  from  measured  peak  pressure  and  impulse. 
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Adjustments  in  the  pressure  load  history  proposed  for  calculating 
response  of  the  ESKIMO  II  east  and  northeast  Igloos  may  be  made  if  more  reliable 
estimates  of  posit Ive-phase  duration  are  obtained  or  if  more  accurate  estimates 
of  the  impulse  are  obtained. 

2.3  PRESSURE  LOAOING  HISTORY  FOR  ESKIMO  IV  TEST 

Eight  pressure  gages  were  placed,  as  shown  in  Figure  2-4,  on  the 
headwall  of  Igloo  B (northeast  igloo)  in  ESKIMO  IV.  This  concentration  of 
instrumentation  was  intended  to  measure  the  variation  of  pressure  over  the 
height  of  the  headwall  and  to  determine  the  uniformity  of  the  pressure  across 
the  face.  In  addition,  side-on  pressure  at  ground  levrl  was  measured  at 
three  ranges  to  establish  the  free-fleld  blast  pressure  corresponding  to  the 
reflected  overpressure  loading  on  the  headwall.  Review  of  the  "qu i ck- look" 
data  indicated  that  useful  records  were  obtained  from  seven  of  the  eight 
stations  on  the  headwall  and  from  all  the  free-field  ground  levei  stations. 

At  the  147-ft  range,  the  free-field  ground  level  station  recorded 
a peak  pressure  of  62  psi.  Empirical  data  for  high-explosive  detonations  in 
air  indicate  that  the  peak  reflected  overpressure  at  the  same  range  should 
be  290  psi  (Ref.  7),  whereas  classical  estimates  of  reflected  shock  strength 
for  ideal  gases  indicate  a peak  reflected  pressure  of  265  psi.  Raw  reflected 
pressure  peaks  taken  from  the  "qui ck- look"  records  range  from  195  to  295  psi. 

Pressure  histories  from  Stations  P1B,  P4B,  and  P5B  appeared  to  be 
influenced  by  movement  or  leakage  of  the  door.  The  record  from  Station  P4B 
suggests  failure  of  the  gage  about  5 msec  after  blast  arrival.  The  record 
from  Station  P78  shows  a prominent  oscillatory  signal  of  about  a 3'msec  period, 
which  distorts  the  early  portion  of  the  record.  Development  of  pressure 
loading  for  the  headwall  of  Igloo  B (northeast)  was  therefore  based  primarily 
on  the  records  from  Stations  P2B,  P3B,  P6B,  and  P8B,  and  on  the  classical 
treatment  and  empirical  data  for  chemical  explosions  and  shock  waves  in  air. 
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FIGURE  2-k.  PLACEMENT  OF  PRESSURE  GAGES,  IGLOO  B (NORTHEAST  IGLOO)  OF  ESKIMO  IV 
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Comparison  of  traces  from  Stations  P2B,  P6B,  and  P8B,  as  shown  in 
Figure  2-5,  fails  to  support  the  assumption  of  systematic  variation  of  the 
pressure  History  over  the  height  of  the  headwall.  Duration  of  the  meaningful 
pressure  loading  on  the  headwall  appears  to  be  14  msec.  This  is  the  duration 
adopted  for  pulses  to  be  used  in  computing  headwall  response. 

Precursors  are  apparent  in  all  three  pulses.  The  delay  between 
earliest  arrival  and  what  is  considered  the  main  blast-pressure  signal  is 
approximately  1 msec.  Ccmyarison  of  signals  from  Stations  P3B  and  P6B  is 
sho*n  in  Figure  2-6.  Although  differences  in  peak  pressure  and  precursor 
timing  and  amplitude  are  apparent,  the  signals  show  good  agreement.  Pressure 
differences  during  the  decay  phase  are  of  the  same  order  as  the  differences 
between  signals  from  the  lower,  middle,  and  upper  stations  shown  in  Figure  2-5* 

On  the  basis  of  the  comparisons  of  "quick-look"  data  for  Stations  P2B, 
P3B,  P6B,  and  P8B,  and  the  analytical  and  empirical  data  for  chemical  explo- 
sions in  air,  a single  pressure  history  was  developed  for  the  headwall  response 
calculation.  The  pulse  is  shown  in  Figure  2-7  with  a smoothed  version  of  the 
free-field  ground  level  (side-on)  pressure  history  measured  at  the  same  range. 
The  peak  pressure,  280  psi,  is  consistent  with  the  analytical,  empirical,  and 
experimental  data.  The  duration,  14  msec,  is  consistent  with  the  ESKIMO  IV 
data,  but  is  considerably  shorter  than  that  predicted  from  the  analytical  and 
empirical  data.  The  impulse,  1180  psi  msec,  is  consistent  with  both  the 
ESKIMO  IV  data  as  shown  in  Figure  2-8  and  with  empirical  data.  The  scaled 
empirical  data  (Ref.  7)  indicate  that  the  reflected  pressure  pulse  correspond- 
ing to  a 62-psi  incident  overpressure  pulse  from  a 37,000-lb  hemispherical 
TNT  charge  should  have  a peak  pressure  of  290  psi,  a pos i t ■ ve-rhase  duration 
of  34  msec,  and  an  impulse  of  1 1 60  psi  msec.  Pulses  closely  resembling  the 
representative  reflected  pulse  shown  in  figure  2-7  can  be  made  to  satisfy 
these  requirements  by  adding  a low  amplitude  tail.  For  example,  the  pulse 
shown  in  Figure  2-9  has  a peak  pressure  of  290  psi,  a duration  of  31*7  msec, 
and  an  impulse  of  1170  psi  msec.  Its  close  resemblance  to  the  proposed  short- 
duration  pulse  is  obvious.  The  use  of  the  short-duration  pulse  is  considered 
satisfactory  and  may  be  preferred  for  reasons  of  economy  in  calculating 
headwall  response. 
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FIGURE  2-6.  COMPARISON  OF  PRESSURE  HISTORIES  TO  TEST  UNIFORMITY  OF  PRESSURE 
LOAOING  ACROSS  UIOTH  OF  WEADWALL,  IGLOO  B OF  ESKIMO  IV 
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FIGURE  2-1.  FREE-FIELD  GROUNO  LEVEL  AHO  REFLECTED  HEADWALL  PRESSURE  PULSES  FOR 
OF  IGLOO  B,  ESKIMO  IV 
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FIGURE  2-8.  COMPARISON  OF  PRESSURE  HISTORIES  FOR  A VERTICAL  ARRAY  OF  GAGES  WITH  THE 
PROPOSED  LOADING  PULSE  FOR  RESPONSE  ANALYSIS  OF  IGLOO  B,  ESKIMO  IV 
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FIGURE  2-9.  COMPARISON  OF  IIEADWALL  PRESSURE  PULSE  FOR  ANALYSIS  OF  IGLOO  6 (ESKIMO  IV) 
WITH  ESTIMATED  REFLECTED  PULSE  ANO  MEASURED  FREE-FIELD  GROUND  LEVEL  PULSE 
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STRUCTURAL  MODEL  OF  HEADWALL 
AND  DOOR  SYSTEM 

This  section  describes  the  finite  element  models  of  the  headwalls 
in  ESKIMO  tests  to  be  used  in  performing  the  dynamic  analyses  described  in 
Section  5*  The  models  are  to  be  used  in  conjunction  with  the  INSLAB  computer 
program.  Section  3.1  describes  the  finite  element  model  of  the  south  headwall 
in  the  ESKIMO  I test.  Also  included  in  this  section  are  the  modeling  tech- 
niques for  considering  the  effects  of  soil  supporting  the  headwall,  concrete 
floor  slab,  and  steel  arch  on  the  response  of  the  headwall.  Section  3.2 
presents  the  finite  element  models  for  the  east  and  northeast  headwalls  in  the 
ESKIMO  II  test.  Section  3*3  discusses  the  finite  element  of  the  northeast 
headwall  in  the  ESKIMO  IV  test. 

3. 1 STRUCTURAL  MODEL.  ESKIMO  I TEST 

3.1.1  ; INITE  ELEMENT  MESH 

Figv  3”1  shows  a coarse  finite  element  mesh  of  the  south  headwall 
in  "•h(?  ESKIMO  I test  which  was  used  in  the  previous  study  to  compute  analyti- 
cal ly  the  dynamic  response  of  the  headwall  (Ref.  2).  A refined  finite  element 
model  of  the  same  headwall  used  in  the  present  study  is  shown  in  Figure  3-2. 

The  previous  model  was  necessarily  coarse  because  the  existing  INSLAB  code 
(Ref.  8)  placed  severe  limits  on  the  size  of  the  model.  To  refine  the  model, 
the  INSLAB  code  was  modified  to  solve  larger  problems  (Ref.  2).  Table  3-1 
presents  a com?.,:  Ison  . ‘he  coarse  and  refined  finite  element  meshes  of  the 
south  headwall  in  ESKIMO  I. 

The  face  of  the  magazine  consists  of  a concrete  headwall  and  two 
steel  wingwalls.  As  di*  ;ed  in  Reference  2,  the  finite  element  model  did  not 
include  the  wingwalls,  ose  effect  on  the  headwalls  is  considered  negligible. 


Preceding  page  blank 
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FIGURE  3-2.  REFI NEC  FINITE  ELEMENT  MODEL  OF  SOUTH  HEABWALL,  ESKIMO  I 


33 


A 


R-7556-1 -4182 


TABLE  3-1 . COMPARISON  OF  COARSE  AND  REFINED 
FINITE  ELEMENT  MESHES 


Item  to  Be  Compered 

Coarse  Mash  Used 
In  Previous  Study 

Refined  Model  of 
Present  Study 

Change, 

% 

1 . 

Total  nodes 

66 

120 

+ 80 

2. 

Total  elements 

51 

96 

+ 90 

3- 

Half  bandwidth  of 
stiffness  matrix 

33 

33 

0 

4. 

Total  elements 
representing 
the  door 

8 

18 

+125 

5. 

• 

Total  discrete 
elements  representing 
soi  1 

43 

71 

+ 65 

6. 

Total  discrete 
elements  representing 
concrete  floor 

0 

* 

18 

-- 

7. 

Nodes  for  which 
response  histories  to 
be  obtained 

10 

19 

+ 90 

8. 

Estimated  fundamental 
frequency  of  the 
largest  element 

1000  Hz 

2250  Hz 

+125 

3** 
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3. 1.1.1  Properties  of  Door 

The  south-igloo  headwal I in  ESKIMO  t was  fitted  with  double-leaf 
hinged  steel  doors.  The  door  was  constructed  from  a 3/o-in.  plate  supported 
by  3“in.  channels  and  3/8-in.  plate  stiffeners,  as  shown  in  Figure  3*3-  A 
No.  16-gage  steel  plate  was  spot-welded  on  the  inside  of  the  door.  The  door 
is  modeled  as  an  orthotropic  plate  of  uniform  thickness  possessing  the  charac- 
teristics of  the  real  door.  The  equivalent  properties  of  the  steel  door  as 
computed  in  Reference  2 are  summarized  below: 


1 


Property 

Symbol 

Numerical  Value 

Equivalent  plate  thickness 

h 

2.04  in. 

Modulus  of  elasticity 

Ex 

32.4  x 106  psi 

Modulus  of  elasticity 

Ev 

27.8  x 106  psi 

Equivalent  density 

p 

0.0812  lb/in.3 

Yield  moment 

M 

XX 

10.52  kip-in. /in. 

Yield  moment 

M 

yy 

12.83  kip-in. /in. 

Plastic  modulus 

ep 

0.1  of  elastic  modulus 
• 

3-1. 2. 2 Properties  of  Hcadwall  Section 

A typical  section  of  the  headwall,  showing  the  details  of  the  rein- 
forcement, is  shown  in  Figure  3”4.  To  account  for  the  rapid  application  of 
loading,  the  dynamic  yield  stress,  f , for  the  reinforcing  steel  is  taken 
as  52  ksi,  as  suggested  in  Reference  Similarly,  the  dynamic  strength  of 
concrete  in  compression,  f^,  is  taken  as  3900  psi,  a 30%  increase  from  the 
static  strength  of  3000  psi.  Young's  moduli  for  steel  and  concrete  are 
assumed  to  be  30  x 10^  psi  and  3 x 10^  psi,  respectively. 
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FIGURE  3-3.  DOOR  DETAILS  FOR  SOUTH  IGLOO,  ESKIMO  I 
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Based  on  Figure  and  the  assumed  material  properties  of 

reinforced  concrete,  the  yield  moments  for  ooth  cross  sections  S and  S 

x y 

are  calculated.  It  is  found  that  both  cross  sections  are  underreinforced; 
therefore,  the  yield  moments  are  governed  by  the  tensile  strength  of  steel. 
The  procedures  for  computing  these  moments  are  presented  in  Appendix  B. 


For  Section 


yy 


121.2  kip-in. /ft  when  the  headwal 1 is  concave  inward 
159.0  kip-in. /ft  when  the  headwal 1 is  concave  outward 


For  Section 


xx 


M 


xx 


110.7  kip-in. /ft  when  the  headwal 1 is  concave  inward 
69.2  kip-in. /ft  when  the  headwal 1 is  concave  outward 


The  headwall  section  is  modeled  as  a homogeneous,  orthotropic  plate 
of  uniform  thickness  possessing  the  same  dynamic  characteristics  as  the  actual 
reinforced  concrete  headwall.  The  equivalent  properties  of  the  headwall 
section  are  given  below: 


Property 

Symbol 

Numerical  Value 

• 

Equivalent  plate  thickness 

h 

12  m. 

Modulus  of  elasticity 

E 

X 

3 x 105  psi 

Modulus  of  elasticity 

E 

y 

2 x 105  psi 

Equivalent  density 

D 

150  lb/ft3 

The  values  of  E and  E are  obtained  usina  the  procedure  described  in 
x y 

Appendix  C. 
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3.'  2 STIFFNESS  OF  SOIL 

The  ESKIMO  tests  (Refs.  1,  4,  5)  showed  that  the  compacted  soils 
behind  the  headwail  and  covering  the  steel  arch  influence  the  behavior  of  the 
headwa 1 1 significantly  because  of  the  interaction  of  soil  and  headwail.  In 
the  longitudinal  direction,  the  earth  cover  merges  smoothly  with  the  ground 
behind  the  rear  wall.  Therefore,  any  support  provided  by  the  soil  per  unit 
area  of  the  headwail  is  similar  to  an  infinitely  long  compressible  column. 

The  response  of  the  soil  is  expected  to  be  highly  nonlinear  because  the 
permanent  displacements  of  the  headwail  are  of  the  order  of  several  inches 
(Refs.  1,  4,  5). 

Actual  test  data  regarding  the  strength  of  the  compacted  soil  over 
the  north  igloo  of  the  ESKIMO  I test  array  are  available  (Ref.  10).  In  the 
test,  two  opposing  rams  were  driven  by  a hydraulic  jack  horizontally  outward 
against  the  earth  fill  through  circular  openings  cut  in  the  right  and  left 
sides  of  the  corrugated  steel  arch.  Load  deflection  curves  were  plotted  for 
the  soil  over  the  igloo  at  the  initial  stage  of  the  loading  and  the  final 
stage  (5  min  after  the  application  of  loads).  The  spring  constants  of  the 
soil  based  on  these  curves  are  shown  in  Table  3*2.  It  is  noted  that  in  the 
testing,  the  soil  was  allowed  to  settle  for  one  minute  before  additional 
loading  was  applied.  For  the  quantities  of  explosives  normally  stored  in 
magazines,  blast-pulse  durations  at  distances  of  interest  are  of  tne  order 
of  tens  of  milliseconds.  As  a consequence,  the  corresponding  stiffness  is 
very  hiah.  To  arrive  at  a reasonable  dynamic  stiffness  of  the  soil,  a 
column  of  soil  subjected  to  step  loading  of  magnitude  P was  invest i gated . 
The  deflection  5 at  time  t is 

PL 

' " AE 
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TABLE  3*2-  SOIL  STIFFNESS  FROM  TEST  IN  NORTH  IGLOO 


K,  lb/in.^ 

Stage 

Locat i on 

P • 0 to  60  psi 

P - 60  to  270  psi 

P > 270  psi 

Right 

700 

350 

230 

Initial 

Left 

666 

245 

133 

Right 

600 

330 

190 

Final 

Left 

530 

235 

133 

The  stiffness  is  defined  as 


|/  _ 

P 

AE 

Iv  • 

6 “ 

L 

where 

A * 

The  cross-sectional  area  of  the 

soil  column 

E - 

PC2 

C - 

Wave  velocity  in  the  soil 

medium,  4000  ft/sec 

p * 

Mass  density  of  the  soil. 

0.00018  lb-sec2/in.4 

L « 

Ct 

Because  of  the  reflections  of  the  pressure  wave  at  the  ends  of  the 
soil  column,  different  expressions  for  the  deflections  were  derived  and  are 
illustrated  in  Figures  3-5*  to  3“5c.  The  resulting  stiffness  with  resoect  to 
time  is  shown  in  Figure  3"5d. 

The  dynamic  stiffness  calculations  are  based  on  a one-dimens ional 
wave  propagation  model.  Neither  soil  friction  nor  wave  dispers;on  was  con- 
sidered; therefore,  the  true  values  should  be  a little  higher  than  the  test  data. 
In  addition,  there  was  ample  time  during  the  test  for  soil  to  settle;  whereas 
in  the  dynamic  stiffness  calculation,  there  is  very  little  tine  for  soil  to 


AO 
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settle.  The  modulus  of  elasticity  for  soil  under  high  pressure  conditions 
is  much  smaller  tnan  that  under  low  pressure  conditions,  as  can  be  seen  in 
Table  3~2,  which  shows  lower  static  stiffnesses  at  high  pressure  levels. 

It  is  doubtful  that  the  effect  of  loading  magnitude  is  important,  especially 
when  the  peak  pressure  occurs  for  only  a 3-msec  duration. 

The  lowest  dynamic  stiffness,  occurring  at  10  msec,  30  msec, 

50  msec,  etc.,  is  864  lb/in.^  (Fig.  3“5d) , which  is  slightly  higner  than  the 
static  stiffness  shown  in  Table  3"2  (700  and  666  lb/in.^  for  right  and  left 
locations,  respectively)  in  the  initial  stage.  Since  the  major  explosive 
loadings  used  in  this  study  were  applied  to  the  headwall  within  the  first 
few  milliseconds,  the  stiffness  is  subject  to  rapid  change  during  this 
period.  A value  of  1300  lb/in. ^ is  considered  to  be  the  best  estimate  on 
the  soil  st i ffness . 

3.1.3  MODEL  REPRESENTATION  FOR  SOIL  AND  HEADWALL  INTERACTION 

In  most  analyses  involving  structures  supported  on  soil,  a common 
practice  is  to  represent  the  soil  by  spring-mass  or  spr i ng-dashpot  systems 
(Ref.  11).  This  approach  is  valid  if  the  motions  are  small  so  that  the  soil 
response  is  linearly  elastic.  Since  the  soil  behind  the  headwall  behaves  ii? 
a highly  nonlinear  fashion,  the  conventional  approach  of  representing  soil 
by  spring-mass  or  spring-dashpot  systems  may  not  be  useful.  In  addition, 
the  presence  of  springs  would  not  allow  the  headwall  to  undergo  permanent 
deflection.  For  this  reason,  viscous  damping  elements,  or  dashpots , were 
used  in  the  previous  study  (Ref.  2).  Massless  dashpots  with  damping 
coefficient  equal  to  cC  per  unit  area  of  headwall  were  used,  where  „ is 
the  mass  density  and  C is  the  wave  velocity  of  the  soil.  The  computed 
displacements  from  the  study  generally  exceeded  the  measured  values  from  the 
test.  The  use  of  dashpots  alone  was  considered  to  be  a factor  in  producing 
excessive  displacements. 
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In  order  tD  determine  the  best  possible  model  to  represent  the  soil, 
the  three  test  cases  shown  in  Figure  3-6  were  considered.  A section  of  the 
headwall  was  sliced  along  the  vertical  plane.  The  base  of  the  headvall  was 
embedded  in  the  ground.  The  top  of  the  arch  line  varies  from  13  ft  at  the 
center  to  zero  at  the  ends  of  the  headwall.  A value  of  5 ft  was  arbitrarily 
select*  ".se  rear  of  the  wall,  from  the  top  of  the  steel  arch  to  the  top 
of  the  i^jadwall,  was  supported  by  earth  fill. 

In  Case  A,  the  earth  fill  was  represented  by  a series  cf  elastic 
springs  attached  to  the  nodes  on  the  headwall,  as  shown  in  Figure  3-6. 

A certain  percentage  of  the  mass  of  soil  was  concentrated  at  each  of  these 
nodes,  as  discussed  in  Reference  2.  The  stiffness  of  the  springs  was  equal  to 
1300  Ib/in.  per  unit  area  of  the  headwall  section.  In  Case  B,  the  earth 
fill  was  modeled  by  a series  of  linear  dampers  attached  to  the  nodes  on  the 
headwall  section.  The  coefficient  of  the  dampers  was  taken  equal  to  2C. 

In  Case  C,  the  soil  medium  behind  the  headwall  was  Tiodeled  by  linear 
dampers  and  nonlinear  springs  in  series  attached  to  the  nodes  on  the  headwall 
section.  The  coefficients  of  the  linear  dampers  were  equal  to  oC.  The 
nonlinear  springs  were  one-way  springs  that  can  resist  notion  only  in  one 
direction,  i.e.,  in  compression.  The  stiffness  of  the  springs  in  compression 
was  equal  to  1300  lb/in.  per  unit  area  of  the  headwall  section. 

All  three  models  were  subjected  to  the  sane  blast-pressure  loading 
experienced  by  the  south  igloo  in  ESKIMO  I.  Dynanic  responses  of  the  models 
in  three  cases  were  obtained  using  TRI/SAC  code  (Ref.  12),  a general  finite- 
element  computer  program  capable  of  solving  linear  and  nonlinear  problems. 

Cases  A and  B were  solved  with  the  linear  step-by-step  option  of  TRI/SAC 
(Ref.  12),  whereas  Case  C was  solved  with  the  nonlinear  step-by-step  option. 

The  response  time  histories  at  selected  locations  are  shown  in 
Figures  3~7  through  3"11-  In  Case  A,  displacement  responses  oscillate  at  high 
frequencies  and  eventually  attenuate  to  zero  (Fig.  3“7).  This  type  of  behavior 
of  the  headwa I 1 is  in  contrast  to  the  observed  behavior  in  the  ESKIMO  I test. 
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FIGURE  3-8.  CASE  B:  DISPLACEMENT-TIME  HISTORIES  OF  THE  HEADWALl 

(Displacements  are  in  inches.) 
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where  lerge  permanent  displacements  of  haedwatls  were  noted.  The  heedwell  in 
Case  S (Fig.  3*8)  does  not  exhibit  any  oscillations  initially  and  then 
stabilizes  at  a displacement  level  slightly  smaller  than  the  maximum  value. 

The  displacement  and  velocity  time  histories  for  Case  C are  shown  in 
Figures  3*10  and  3“11»  respectively.  The  behavior  of  the  headwall  in 
Case  C is  similar  to  Case  B except  that  the  final  displacements  in  Case  C 
are  well  below  the  maximum  values. 

In  the  ESKIMO  I test,  a gap  of  several  inches  between  the  headwall 
and  the  earth  fill  was  noted  (Ref.  1).  This  indicated  that  the  permanent 
displacements  of  the  headwalls  were  significantly  smaller  than  the  maximum 
values  reached  during  the  test.  The  soil  model  of  Case  C (Fig.  3~&)  appears 
to  best  represent  the  qualitative  behavior  of  the  earth  fill  in  the 
ESKIMO  tests.  However,  the  magnitude  of  the  displacements  obtained  in 
Case  C indicates  that  reduced  numerical  values  of  coefficients  for  the  springs 
and  dampers  must  be  used.  In  the  absence  of  any  relevant  measurements,  the 
only  guidelines  for  choosing  these  coefficients  is  their  performance  in  small 
test  problems.  It  was  found  that,  in  such  a test  problem,  a spring  constant 
of  40  lb/in.  per  unit  area  and  a damping  coefficient  of  0.75  lb-sec/in.  per 
unit  area  produced  headwall  displacements  that  were  of  the  order  of  magnitude 
observed  in  ESKIMO  I. 

This  combination  of  linear  damper  and  nonlinear  one-way  spring  was 
used  to  represent  the  interaction  of  soil  with  the  finite  element  mesh  of 
the  south  headwall  (Fig.  3*2).  The  constants  for  the  spring  and  dashpot  at 
a node  were  based  on  the  total  tributary  area  of  soil  to  be  concentrated  at 
that  node.  The  existing  INSLAB  code  did  not  provide  for  nonlinear  springs. 
Therefore,  the  program  was  modified  to  accept  nonlinear  one-way  springs. 
Theoretical  details  on  the  incorporation  of  these  springs,  along  with  some 
test  cases,  are  presented  in  Appendix  A. 
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3.1.4  EFFECTS  OF  CONCRETE  FLOOR  SLAB 

The  floor  of  the  south  Igloo  in  ESKIMO  I is  a concrete  slab  6-in. 
thick,  which  thickens  to  12  in.  near  the  headwall.  The  reinforcement  consists 
of  No.  k ba:s  spaced  12  in.  apart,  3 in.  below  the  top  surface  of  the  slab. 
Linear  spring  elements  and  concentrated  masses  are  used  to  model  the  effects 
of  the  floor.  The  spring  constant  per  unit  area  of  the  floor  slab,  determined 
in  response  to  a static-compressive  stress,  is 


f 

4 


kL  - 12,500  lb/in.3 


This  value  was  obtained  assuming  Young's  modulus  of  concrete  to  be 

E - 3 x 106  lb/in.2 

c 


and  the  length  of  the  igloo  to  be  20  ft.  The  stiffness  of  the  floor  in 
response  to  a dynamic  load  would  be  larger  than  the  value  given  above. 
However,  the  fact  that  the  rear  wall  of  the  igloo  is  not  anchored  perfectly 
would  result  in  a slightly  lower  stiffness  value.  Assuming  that  the  various 
unknowns  compensate  for  each  other,  the  static  stiffness  was  used  in  modeling 
the  effect  of  the  floor. 

The  value  of  the  mass  of  the  floor  slab  participating  in  the  mass/ 
spring  representation  was  chosen  to  preserve  the  fundamental  frequency  of 
longitudinal  vibration  of  the  floor.  The  resulting  value  is  equal  to  the 
mass  of  a strip  of  floor  slab  3.2-ft  wide. 
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3.1.5  EFFECTS  OF  STEEL  ARCH 

The  main  function  of  the  steel  arch  in  a magazine  is  to  protect  the 
contents  of  the  magazine  from  ti.e  weather  and  from  the  earth  fill  over  the 
magazine.  The  only  structural  requirement  for  the  arch  is  to  support  the 
earth  fill  around  the  arch.  The  material  used  in  the  arch  for  the  south 
igloo  headwal 1 (ESKIMO  1)  is  No.  1-gage  corrugated  steel.  The  corrugation 
increases  the  bending  and  circumferential  stiffness  of  the  arch,  but  reduces 
its  longitudinal,  or  axial,  stiffness. 

To  study  the  longitudinal  stiffness  of  the  arch,  one  corrugation  of 
the  arch  is  considered.  Since  this  is  symmetrical  at  the  centerline  of  the 
foU  only  one  half  of  the  f^ld  is  necessary.  The  curved  fold  is  approximated 
by  a straight  segment  AB,  as  chown  in  Figure  3~12.  Under  the  action  of  a 
horizontal  load  P at  end  A of  the  segment,  the  horizontal  displacement  of 
A results  due  to  bending  and  axial  compression  of  the  member  AS.  The 
component  of  deflection  due  to  bending  can  be  calculated  as 

5 s PL3  sin2  9 

b 12E I cos3  0 

and  the  component  due  to  axial  compression  is  given  by 

. _ PL  cos  0 

°a  " AE 

where  El  and  AE  are,  respect i vely , the  bending  and  axial  stiffnesses  of 
the  arch.  Therefore  the  total  horizontal  displacement  of  A is 

, _ PL3  sin2  3 PL  cos  0 

12EI  cos 3- 
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For  a unit  width  of  the  arch,  the  following  values  are  used  to 

compute  : 

A - 0.2758  in.2 

E * 30  x 10^  psi 
I - 0.001 7^3  in.4 

L ■ 3 »n. 
h - 0.2758  in. 


The  horizontal  displacement  5 is  computed  as 
6 - (2.32  x 10-5  in. /lb)  P 


If  the  corrugated  arch  were  to  be  treated  as  if  it  were  constructed 
of  an  equivalent  uncorrugated  material,  the  relation  between  load  and  deflec- 
tion would  be 


. PL_ 

2 * AE 

e 

where  is  the  equivalent  elastic  modulus  of  the  uncorrugated  material. 

This  relation  can  be  used  to  determine  the  ec-lvalent  moculus,  making  ^se 
of  the  above  relation  between  load  and  aef lection.  The  result  is 

E * 4.6G  x 10^  psi 

e 

The  buckling  stress  of  tne  equivalent  uncorrugates  arc-  is  cc'c-tec 
form  the  relation 


E h 


e 


5* 
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where  a is  the  radius  of  the  arch  (a  - 156  in.),  h is  the  thickness  of  the 
arch,  and  v is  Poisson's  ratio  of  the  arch  material  (assumed  -j  ■ 0.25) ■ 

The  buckling  stress  is  computed  to  be  495  psi.  Although  this  is  the  theo- 
retical buckling  stress  for  the  equivalent  uncorrugated  arch,  buckling 
experiments  on  cylindrical  shells  show  that  the  stress  at  which  the  shells 
buckle  is  seldom  better  than  30%  of  the  theoretical  value  (Ref.  13).  Thus, 
the  buckling  stress  of  the  arch  is  estimated  to  be  150  psi. 

Nodal  point  31  of  the  finite  element  grid  of  Figure  3*2  is  a 
typical  nodal  point  located  on  the  arch.  The  length  of  the  arch  that  is 
attributed  to  tnis  node  is  about  24  in.  Therefore,  the  force  acting  on  this 
node  due  to  buckling  of  the  arch  is  990  lb.  At  no  time  can  the  force  of  the 
arch  on  Node  31  exceed  this  value.  This  node  is  also  influenced  by  an  area 
of  about  143  in . ^ of  soil  backfill.  Using  Node  28  (see  Fig.  3" 1 ) of  the 
previous  calculation  of  ESKIMO  I,  south  igloo,  as  a guide,  the  expected 
peak  displacement  and  velocity  of  Node  31  are  5 in.  and  150  in. /sec, 
respectively.  The  soil  is  represented  with  springs  having  a stiffness  of 
40  Ib/in.  per  unit  area,  and  dampers  having  a viscosity  of  0.75  )b-sec/in. 
per  unit  area.  These  values  result  in  a peak  force  of  28.600  lb  due  to 
elasticity  of  the  soil,  and  an  additional  peak  force  of  16,1.00  lb  due  to 
damping.  By  comparison,  the  contribution  of  the  arch  is  expected  to  be 
negl igible. 

A study  of  the  permanent  deformations  and  damage  patterns  of  headwalls 
in  the  ESKIMO  tests  indicate  that  the  steel  arch  appears  to  strongly  influence 
the  behavior  of  the  headwalls  during  the  tests.  As  an  example,  Figure  3*13 
is  a contour  diagram  of  the  deformation  of  a headwall  after  the  ESKIMO  I 
test.  The  deformation  contours  generally  approximate  the  shape  of  the  steel 
arch,  which  was  semicircular.  Figure  3” 1 4 shows  a front  photographic  view 
of  the  same  headwall  after  the  test.  The  crack  pattern  experienced  by  the 
headwall  again  approximates  the  shape  of  the  steel  arch.  It  is  believed 
that  this  phenomenon  occurs  because  the  arch  separates  the  supported  and 
unsupported  areas  of  the  headwall.  In  particular,  the  unsupported  area  of 
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FIGURE  3 13-  MOVEMENT  OF  HEAPWALL  OF  NORTH  ACCEPTOR  IGLOO,  ESKIMO  I. 
A plus  value  shows  movement  away  frvm  the  donor 
magazine;  a minus  value  shows  movement  toward. 

The  units  are  in  hundredths  of  feet. 
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the  headwell  is  semici rcular,  and  hence  it  is  not  surprising  that  the  permanent 
deformation  and  damage  patterns  of  the  headwell  resemble  the  shape  of  the  arch. 
In  addition,  the  sharpness  of  the  steel  arch  causes  high  stresses  in  the  head- 
wall,  causing  local  failure  along  the  arch  line. 

3-2  STRUCTURAL  MODELS  OF  HEADWALLS  IN  ESKIMO  II 

Two  different  headwalls  in  the  ESKIMO  I I test  have  been  analyzed. 
These  are  the  Igloo  C (east)  and  Igloo  B (northeast)  headwalls.  The  structural 
models  of  these  headwalls  are  described  in  the  following  subsections. 

3.2.1  IGLOO  C (EAST)  HEAOWALL 

As  described  in  Section  1.2,  the  east  headwall  in  the  ESKIMO  II  test 
is  identical  in  every  respect  to  the  igloos  used  in  ESKIMO  I.  Therefore,  the 
finite  element  model  of  the  east  headwall  in  ESKIMO  II  is  identical  to  the 
model  of  the  south  headwall  in  ESKIMO  I (Fig.  3“2)  described  in  Section  3-1. 

3.2.2  IGLOO  B (NORTHEAST)  HEADWALL 

Figure  3-15  shows  plan  and  an  elevation  of  the  Igloo  B (northeast) 
headwall  in  the  ESKIMO  II  test.  The  wingwalls  are  of  reinforced  concrete 
construction  and  are  monolithic  with  the  headwall.  Therefore,  the  finite 
element  model  of  the  headwall  shown  in  Figure  3"l6  includes  the  portion  of 
the  wingwalt  up  to  the  counterfort  in  order  to  consider  the  effect  of  the 
wingwall  on  the  response  of  the  headwall.  As  described  before,  an  oval-shaped 
steel  arch  is  used  behind  the  headwall  to  support  the  soil.  The  opening  in 
the  headwall  for  the  door  is  strengthened  with  concrete  beams,  as  shown  in 
Figure  3“ 15-  On  top  of  the  horizontal  beam,  there  is  an  overhang  that  extends 
3 ft  from  the  headwall.  Since  there  is  very  little  reinforcement  provided 
along  the  vertical  direction,  the  overhang  is  not  considered  to  function  as  a 
structural  member.  Therefore,  although  the  mass  of  the  overhang  is  included 
in  the  finite  element  model,  the  stiffness  is  not.  In  designing  the  finite 
element  mesh  shown  in  Figure  3'l6,  the  beams  are  represented  by  separate 
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FIGURE  3-15.  IGLOO  C (NORTHEAST)  HEAOMALL.  ESKIMO  II 


FIGURE  3-15.  (CONCLUDED) 


elements  so  that  the  actual  beam  properties  can  be  assigned  to  only  the 
elements  representing  the  beams.  This  avoids  the  averaging  of  properties 
that  becomes  necessary  when  the  same  element  represents  the  beam  and  a 
portion  of  headwal 1 . The  arch  line  is  placed  as  closely  as  possible  on  the 
edges  of  the  finite  elements.  In  order  to  minimize  the  use  of  triangular 
elements,  the  arch  line  is  permitted  to  run  across  the  Elements  11,  12,  22, 

33*  and  44,  as  shown  in  Figure  3“16.  The  finite  element  shown  in  Figure  3-1& 
consists  of  a total  of  107  plate  elements  and  127  nodes. 

1 The  biparting  and  sliding  steel  door  of  the  northeast  igloo  consists 

of  a 5/8~in.  plate  and  a 1/4-in.  plate  sandwiching  five  51 10  steel  sections 
arranged  vertically  as  shown  In  Figure  3“17.  Following  the  procedures 
described  in  Reference  2,  the  material  properties  of  the  elements  representing 
the  steel  door  are  computed.  These  values  are  summarized  below: 


Property 

Symbol 

Numerical  Value 

Equivalent  plate  thickness 

h 

4.14  in. 

Modulus  of  elasticity 

E 

30  x 103  ksi 

Equivalent  densi ty 

p 

0.0773  lb/in.3 

Yield  moment 

M 

XX 

68.6  kip-in. /in. 

Yield  moment 

M 

yy 

92.3  ki p- in . /i n . 

Plastic  modulus 

E 

P 

0.1  of  elastic  modulus 

Poisson's  ratio 

V 

0.25 

The  headwal 1 and  wingwall  are  modeled  as  homogeneous  plates  of 
uniform  thickness  possessing  the  same  dynamic  characteristics  as  the  actual 
rei nforced-concrete  sections.  The  yield  moments  for  these  plates  are  computed 
using  the  procedures  described  in  Appendix  B.  The  yield  moments  and  equivalent 
properties  of  the  headwal 1 and  wingwall  sections  are  shown  in  Table  3~3- 
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TABLE  3-3.  PROPERTIES  OF  HEAOWALL  AND  WINGWALL  SECTIONS 
OF  IGLOO  B (NORTHEAST),  ESKIMO  II 


Property 

Value  for 
Headwall  Section 

Value  for 
Wingwall  Section 

Equivalent  plate  thickness,  h 

12  In. 

12  in. 

Modulus  of  elasticity,  E 

x 1U5  psi 

2*67  x 10^  psi 

Equivalent  density,  o 

ISO  lb/ft3 

ISO  lb/ft3 

Yield  moment,  M^ 

290  kip-tn./ft 

101  kip-in. /ft 

Yield  moment,  M^ 

-231  kip-in. /ft 

-79  kip-in. /ft 

Yield  moment,  M*y 

196  kip-in. /ft 

192  kip-in. /ft 

Yield  moment,  M~y 

-160  kip-in. /ft 

-83  kip-in. /ft 

AA8012 


The  supporting  system  of  the  headwall  consists  of  the  backfill  on 
top  of  the  steel  arch,  the  concrete  floor  slab,  and  the  soil  below  the  fl cor 
slab.  These  supports  are  represented  in  the  form  of  springs  and  dampers  as 
described  in  Sections  3-1-3  and  3-1-**- 


3-3  STRUCTURAL  MODEL  OF  IGLOO  B (NORTHEAST)  HEADWALL,  ESKIMO  IV 

Figure  3-18  shows  the  finite  element  model  of  the  headwall  of  tre 
northeast  igloo  in  the  ESKIMO  IV  test.  An  oval-shaped  steel  arch  suooorts 
the  soil  for  this  igloo.  The  door  opening  in  the  headwall  is  s ■;  re-c  tne-eo 
with  concrete  beans.  The  vertical  Dean  extends  to  the  too  of  the  neaa..al’. 

As  in  the  case  of  the  northeast  headwall  in  tne  ESKIMO  II  test,  t-e  -cdei  of 
this  headwall  represented  the  beams  by  separate  elements.  An  e*fcrt  was  ~ade , 
as  far  ar  possiole,  to  place  tne  arch  line  on  the  edges  cf  t-e  - i - i t e e'e-e-ts 
However,  in  order  to  mini-.ice  the  use  of  triangular  eie“ent5.  t-e  a-;-  !'-e 
was  permitted  to  run  across  Elements  22.  33.  **3.  ana  53.  as  sncwr  in 
Figure  3“ 1 S . The  finite  element  mes*  ior  this  head-all  s c-ec  so  rat 


A 


S 10 


nnnouiiiaii 


73  74  75 


•2  93  94 


91  92  93 


100  101  102  103 


109  110  111  112 


BO.  81 


89  /90 


98  / 99 


107  / 108 


116  117 


118  119  UQ  121 


is  h ■ 


125  126 


FIGURE  3- '3.  FINITE  ELEMENT  MODEL  OF  NORTHEAST  HEADWALL.  ESKIMO  IV 
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the  locations  of  the  accelerometers  and  the  LVDT  transducers  (see  Sec.  k)  gen- 
erally correspond  to  seme  nodes  in  the  mesh.  This  facilitates  a direct  compari- 
son of  the  test  and  calculated  results.  The  finite  element  mesh  of  the  head- 
wall  shown  in  Figure  3“l8  consists  of  a total  of  126  elements  and  150  nodes. 

The  single- leaf  sliding  steel  door  of  this  headwall  consists  of  a 
3/8-in.  plate  and  a 1/4-in.  plate,  strengthened  by  four  6WF25  wide  flange 
sections  and  two  6 C 13  channels,  all  running  horizontally.  Two  6 £ 13  channels 
also  run  vertically  along  the  two  ends.  This  arrangement  is  shown  in  Figure  3-19- 
Following  the  proceduros  described  in  Reference  2,  the  material  properties  of 
the  elements  representing  the  steel  door  were  computed.  These  values  are 
summarized  below: 


Property 

Symbol 

Numer i ca 1 Va 1 ue  j 

Equivalent  plate  thickness 

h 

k.kS  in. 

Modulus  of  elasticity 

E 

30  x n3  ksi 

Equivalent  density 

0 

0.05;-  lb/in.3 

Yield  momen  t 

M 

XX 

107-0  ki p- i n . / in . 

Yield  moment 

M 

9^.5  ki p- in . / in . 

VY 

Plasti c modul us 

E 

0.1  of  elastic  modulus 

P 

Poisson's  ratio 

V 

0.25 

To  account  for  the  variation  of  section  properties,  a total  of  nine 
different  materials  is  used  to  describe  the  model,  as  shown  in  Figure  3“2Q. 
Material  No.  1 refers  to  the  properties  of  the  steel  door,  as  discussed  above. 
Due  to  the  different  a r,  ang  intent  of  the  reinforcing  steel,  the  headwall  section 
is  assigned  four  different  materials  with  Nos.  2,  3,  8,  and  9,  as  shown  in 
Figure  3“20.  Material  Nos.  5,  6,  and  7 are  used  to  represent  the  beams. 
Material  '<  is  assigned  r~  trie  footing  of  the  headwall.  The  material 
properties  assigned  to  each  material  number  is  shown  in  Table  The  yielc 

moments  of  ail  concrete  sections  are  computed  using  the  procedures  described 
in  Appendix  B. 


(a)  Elevation 


(b)  Longitudinal  section 


(c)  Cross  section 


figure  3-19.  configuration  of  steel  door  for  igloo  b (northeast)  headwall, 

ESKIMO  IV 


FIGURE  3-"'0.  REGIONS  OF  DIFFERENT  MATERIAL  PROPERTIES  FOR 
IGLOO  B (NORTHEAST)  HEADWALL,  ESKIMO  IV 
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The  supporting  system  of  the  headwall  consists  of  tne  backfi 
top  of  the  steel  arcn,  the  concrete  floor  slab,  ana  the  soil  below  the 
slab.  The  representat ion  of  these  supports  in  the  form  of  springs  and 
dampers  is  described  in  Sections  3-1 -3  and  3-1-**. 
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SECTION  4 

INSTRUMENTATION  FOR  ESKIMO  IV 

This  section  describes  Agbabian  Associates'  recommendations  on  the 
i ns t rumen  tat  ion.  of  the  igloos  in  the  ESKIMO  IV  test:  types  of  transducers  fiat 

*»ere  to  be  used,  their  locations  on  the  headwalls  and  doors,  and  tne  ranees  oe 
variables  over  which  the  transducers  were  to  be  calibrated.  The  recommenda- 
tions reflected  a desire  to  obtain  sufficient  measurements  of  the  northeast 
headwall,  since  this  is  the  only  headwall  of  ESKIMO  IV  that  Mas  analyzed 
using  tne  INSLAB  code. 


4.1  PRESSURE  GAGES 

The  Kistler  gages  were  used  to  measnr * pressure- t ime  histories  on 
headwalls  in  the  ESKIMO  IV  test.  These  quartz  crystal  piezoelectric  trans- 
ducers were  selected  because  of  their  excellent  past  performance  in  similar 
events.  Table  4-1  shows  estimated  peak  pressures  and  the  ratings  of 
pressure  gages  that  were  scheduled  to  be  used  in  ESKIMO  IV.  Table  4-1  also 
shows  similar  information  for  the  Ballistic  Research  Laboratory  (BRL)  self- 
recording  mechanical  gages,  placed  in  pairs,  to  measure  far-field  pressures 
along  radial  lines  in  tho  northwest,  west,  and  south  directions. 


t 
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Figure  4-1  shows  the  recommended  arrangement  of  pressure  gages  on 
the  northeast  igloo  (Igloo  B) . The  primary  objective  of  this  recommendation 
was  to  determine  the  variation  of  pressure  loading  on  the  wall  in  vertical  and 
horizontal  directions.  Because  of  the  symmetry  of  the  loading  with  respect  to 
the  centerline  of  the  headwall,  all  gages  except  one  were  placed  on  the  right 
half  of  the  headwall.  The  gage  on  the  left  half  of  the  headwall  was  intended 
to  detect  any  significant  deviations  in  the  symmetry  of  the  loading.  The 
arrangement  shown  in  the  figure  commits  8 pressure  gages  for  the  northeast 
igloo.  Figure  4-2  shows  the  recommended  arrangement  of  the  pressure  gages 
on  the  remaining  two  igloos. 
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4.2  ACCELEROMETERS 

Accelerometers  (Model  114  by  Setra  Systems  Inc.)  were  used  to 
measure  acceleration- time  histories  at  various  locations  on  headwalls  in 
the  ESKIMO  IV  test.  These  gages  yield  good  response  from  zero  hertz  to 
several  hundred  hertz,  and  the  damping  of  the  seismic  element  is  essentially 
independent  of  the  changes  in  ambient  temperature.  Table  4-2  shows  estimated 
peak  accelerations  and  ratings  of  the  accelerometers  used  in  ESKIMO  IV. 

Recommended  placement  of  the  accelerometers  on  the  headwall  of  the 
northeast  (B)  igloo  is  shown  in  Figure  4-3.  Because  of  the  assumed  symmetry  of 
the  loading  and  the  structure,  the  transducers  were  placed  only  on  one  side 
of  the  centerline  of  the  igloo.  The  recommended  locations  for  the  gages 
generally  coincided  with  the  nodes  of  the  finite  element  mesh  that  was 
used  to  analyze  the  headwall.  As  shown  in  the  figure,  a total  of  12  acceler- 
ometers was  placed  on  this  headwall  and  no  accelerometers  were  placed  on  the 
headwalls  of  the  east  (0)  and  west  (E)  igloos. 

The  impulse  response  of  the  accelerometers  was  checked  out  prior  to 
the  test  with  a drop-table  calibration  device.  This  portable  device  was 
carried  to  the  site  for  one-the-spot  checkout  of  the  accelerometers. 

4.3  LINEAR  MOTION  TRANSDUCERS 

Long-stroke  displacement  transducers  of  the  LVDT  type  were  used  in 
ESKIMO  IV  to  measure  movement  of  the  concrete  headwall.  The  female  part  of 
the  transducer  was  attached  to  a universal  joint  mounted  on  a rail  section 
(mass)  that  was  suspended  from  the  corrugated  metal  roof  of  the  igloo  by 
chains.  The  male  or  rod  portion  was  supported  by  universal  joints,  with  one 
universal  joint  securely  attached  to  the  headwall  at  the  desired  point  of 
measurement.  An  effort  was  made  to  reference  as  many  transducers  as  practical 
to  the  same  mass. 
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FIGURE  4-3.  PLACEMENT  OF  ACCELEROMETERS  ON  HEADWALL  OF  NORTHEAST 
IGLOO,  ESKIMO  IV 
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TABLE  k-2.  SCHEDULE  OF  ACCELEROMETERS  FOR 
NORTHEAST  IGLOO.  ESKIMO  IV 


Position 

Estimated 
Maximum 
Accelerat ion , 
9 

Oesi red 
Accelerometer 
Rating, 

9 

Door 

630 

1000 

Headwal)  --  away  from  door 

150 

300 

Headwall  — near  door 

200 

300 

Figure  k-k  ihows  recommended  placement  of  the  displacement  trans- 
ducers (LVDT)  on  the  headwal 1 of  the  northeast  igloo.  Again,  the  locations 
shown  for  LVOTs  generally  coincided  with  the  nodes  of  the  finite  element  mesh. 
Nine  LVOTs  were  placed  on  the  northeast  igloo,  as  shown  in  the  figure,  and 
6 LVOTs  were  placed  or  the  remaining  two  Igloos. 
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FIGURE  k-k.  PLACEMENT  OF  LVOT  TRANSDUCERS  ON  HEADWALl  OF  NORTHEAST 
IGLOO,  ESKIMO  IV 
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SECTION  5 

DYNAMIC  RESPONSE  OF  HEAOUALLS  IN  ESKIMO  TESTS 

5.1  INTRODUCTION 

This  section  describes  the  results  of  dynamic  analyses  of  four 
different  headwalls  in  tne  ESKIMO  tests.  The  four  headwalls  considered  are 
the  south  (and  west)  headwall  in  the  ESKIMO  I test,  the  east  and  northeast 
headwalls  in  the  ESKIMO  II  test,  and  the  northeast  headwall  in  the  ESKIMO  IV 
test.  The  analyses  were  performed  subjecting  the  finite  element  model  of  the 
headwall  and  door  systems  to  blast  loads.  The  finite  element  models  are 
discussed  in  Section  3 and  the  pressure  loadings  applicable  to  each  model 
are  given  in  Section  2.  The  dynamic  response  of  the  models  was  obtained 
with  the  modified  INSLAB  code  (Ref.  3) > a nonlinear,  dynamic  finite  element 
program.  The  results  of  these  analyses  are  compared  with  the  available  test 
results  from  the  ESKIMO  tests. 

Because  Section  5 requires  some  70  figures  to  present  the  analysis 
results,  they  have  been  grouped  at  the  end  of  the  section  text.  Figure  5-1 
et  seq.  begins  on  page  106. 

• 

5.2  RESPONSE  OF  SOUTH  HEADWALL,  ESKIMO  I 

The  results  of  the  ESKIMO  I calculation  using  the  refinid  finite 
element  model  of  the  south  headwall  (Fig.  3"2,  Sec.  3)  are  presented  in 
this  section.  The  finite  element  model  of  the  headwall  was  subjected  to 
the  pressure  loadings  shown  in  Figure  2-1  (Sec.  2),  The  results  of  the 
calculation  are  also  applicable  to  the  west  headwall,  since  the  south  and 
west  headwalls  have  identical  geometries  and  were  exposed  to  the  same  loading. 
These  results  are  shown  as  time-history  plots  of  the  headwall  motion  and 
contour  plots  of  the  displacement  of  the  headwall.  Comparisons  are  made 
between  these  results  and  the  results  of  the  prior  calculation  of  Reference  2. 
In  addition,  where  experimental  measurements  are  available,  comparisons  are 
made  between  the  calculation  and  the  actual  response  of  the  south  and  west 
igloos  of  ESKIMO  I . 
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5.2.1  TIME-HISTORY  DATA 

Figures  5-1  through  5-9  show  displacement,  velocity,  and  accelera- 
tion histories  computed  at  several  locations  on  the  south  and  west  headwalls 
and  doors  of  the  ESKIMO  I event.  In  addition,  Figures  5-10  through  5-12  show 
displacement  histories  at  several  other  locations.  Typically,  these  response 
histories  exhibit  a steady  climb  to  peak  displacement,  followed  by  a recovery 
of  displacement  such  that  the  magnitude  of  the  final  displacement  is  less 
than  half  the  magnitude  of  the  peak  displacement.  Furthermore,  the  frequency 
content  of  the  displacement  histories  is  extremely  low  compared  to  the  input 
pressure  histories.  The  rise  times  of  the  displacement  hismries  are  often 
30  or  40  msec,  whereas  the  rise  time  of  the  input  pressure  history  is  1 msec. 
This  is  because  the  fundamental  frequency  of  the  headwall  is  very  low 
compared  tc  the  predominant  frequencies  in  the  input  loading. 


The  tendency  for  final  displacements  to  be  smaller  than  the  peak 
displacements  is  due  primarily  to  the  use  of  nonlinear  spring  elements  In 
modeling  the  soil  backfill  behind  the  headwall.  A second  factor  is  the  use 
of  linear  springs  to  model  the  concrete  floor  slab.  In  contrast,  the 
calculations  of  Reference  2 used  only  damping  elements  to  model  both  the 
backfill  and  the  foundation.  This  latter  formulation  resulted  in  final 
displacements  that  were  nearly  equcl  to  the  peak  d:splacements. 

Figure  5-1  shows  the  motion  of  Node  1,  which  corresponds  to  the 
top  center  of  the  headwalls  of  the  south  and  west  Igloos  of  ESKIMO  I.  The 
peak  positive  displacement*  is  0.3  in.,  whereas  the  final  displacement  is 
negative  and  is  equal  to  0.7  in.  This  indicates  the  presence  of  a gap  about 


*ln  this  '•epcrt  the  positive  displacement  of  a headwall  represents  movement 
away  from  the  source,  or  into  the  igloo;  and  the  negative  displacement 
represents  movement  towards  the  source,  or  out  of  the  igloo. 
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1-ir.  wide  between  the  headwal]  and  the  compressed  backfill  after  the  shock. 
Figure  5“13  (from  Ref.  2),  which  shows  the  motion  of  the  corresponding  point 
of  the  prior  calculation,  predicts  no  such  gap.  However,  gaps  between  the 
headwal 1 and  the  soil  have  been  observed  during  posttest  inspection  of  the 
acceptor  magazines.  This  indicates  that  the  present  method  of  modeling  the 
soil  is  a step  in  the  right  direction.  The  displacement  histories  of  other 
nodes  on  the  backfill  portions  of  the  headwal Is  are  shown  in  Figure  5*10. 

Figure  5*2  shows  the  response  history  of  Node  19,  located  on  the 
top  of  the  arch  line  of  the  south  and  west  headwal Is.  The  response  of  the 
corresponding  point  of  the  prior  calculation  is  shown  in  Figure  5*1^  (from 
Ref.  2).  The  displacement  histories  begin  in  very  similar  fashion,  reaching 
displacements  of  several  inches  after  1*0  msec.  At  that  time  the  present 
calculation  begins  to  recover  displacement,  whereas  the  prior  calculation 
continues  to  show  increased  displacement.  A gap  of  about  2.75  in.  develops 
in  the  current  calculation.  A large  difference  in  final  position  of  the 
headwall  is  shown  in  the  two  calculations.  The  current  calculation  shows 
a final  displacement  of  about  0.3  in.,  whereas  the  prior  calculation  shows 
a 6-in.  final  displacement.  A similar  situation  exists  with  Figure  5*1. 

Figure  5*1 » shows  the  displacement  histories  of  several  other  locations 
along  the  arch  line. 

Figures  5-5  and  5*8  show  the  responses  of  points  located  within 
the  arch  line  on  the  south  and  west  headwal Is.  While  less  strongly  influenced 
by  the  soil  backfill  than  points  located  on  or  outside  the  arch  line,  these 
points  also  show  considerable  recovery  from  peek  displacement.  However, 
the  ratios  of  final  displacement  to  peak  displacement  are  greater  in  these 
cases  than  in  those  more  directly  dependent  on  soil  behavior.  Figure  5*15 
(from  Ref.  2)  shows  the  response  for  the  prior  calculation  of  the  por"* 
corresponding  to  that  of  Figure  5*8.  The  displacement  histories  bet.,  ,i 
the  two  cases  disagree  in  the  ratio  of  final  displacement  to  peak  displace- 
ment, and  in  the  magnitude  of  the  peak  displacement.  The  two  cases  are  in 
good  agreement  on  the  peak  velocity  and  the  peak  acceleration,  although  the 
durations  are  different. 
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Figures  5-3,  5-4,  and  5*7  show  the  responses  of  points  on  the 
boundary  between  the  headwall  and  the  steel  door.  In  the  present  calculation, 
the  door  Is  connected  to  the  headwall  only  along  its  top  and  side  edges.  It 
is  free  along  its  bottom  and  center  edges.  As  a result,  the  door  exerts  a 
great  deal  of  influence  on  the  headwall  along  the  contact  boundary.  Since 
the  connection  along  this  boundary  is  accomplished  by  hinges,  all  of  this 
influence  is  transmitted  to  the  headwall  in  the  form  of  shear  force.  The 
resulting  displacements  are  considerably  larger,  therefore,  than  displace- 
ments of  the  headwall  not  along  the  boundary.  For  example,  the  displacement 
of  Node  49  (Fig.  5-4)  exceeds  that  of  the  adjacent  point.  Node  50  (Fig.  5-5). 
Similarly,  the  displacement  of  Node  76  (Fig.  5-7)  exceeds  that  of  Node  77 
(Fig.  5-8).  Figures  5-16,  5-17,  and  5-18  (from  Ref.  2)  show  points  of  the 
prior  calculation  that  correspond,  respectively,  to  Figures  5-3,  5-4,  and 
5-7  of  the  current  calculation.  In  the  previous  calculation,  the  door  was 
modeled  as  having  shear  support  along  the  bottom  boundary  as  well  as  along 
the  top  and  side.  As  a result,  the  door  played  a slightly  less  important 
role  in  the  displacements  of  the  prior  calculation.  The  peak  displacements 
of  the  present  calculation  are  comparable  to  those  of  the  prior  calculation; 
and,  in  the  case  of  Node  46  (F:g.  5“3) , the  peak  displacement  greatly 
exceeds  that  of  the  corresponding  point  in  the  prior  calculation.  However, 
the  final  displacements  are  below  those  of  the  prior  calculation. 


Experimental  data  from  the  south  igloo,  for  a location  corresponding 
to  Node  46,  is  available.  A comparison  between  the  measured  data  and  the 
computed  response  is  given  in  Table  5-1.  Although  the  average  acceleration 
is  in  good  agreement,  the  duration  and  the  velocities  are  not.  The  trend  of 
the  comparison  suggests  that,  locally,  the  door  may  have  broken  free  of  its 
support  during  loading.  Had  this  been  able  to  occur  in  the  finite  element 
model,  the  unloading  of  Node  46  would  have  been  initiated  much  earlier, 
resulting  in  lower  velocities.  This  would  have  been  accomplished  without 
significantly  altering  the  sverage  acceleration.  It  would  be  desirable  in 
future  calculations  to  possess  the  capability  of  modeling  breakable  hinges. 
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TABLE  5-1.  MOTION  OF  SOUTH  HEADWALL  (ESKIMO  I) 
ABOVE  CENTER  OF  DOOR 


1 tern 

Measured  Value 
for  South  Igloo 

Present  Calc. 
(Node  46) 

Prior  Calc. 
(Node  33) 

Maximum  Velocity,  fps 

27-9 

44.2 

20 

Average  Velocity  from 
Initial  Motion  to  Peak 
Displacement,  fps 

13-3 

29.0 

11.7 

Time  Interval  from 
Initial  Motion  to 
Peak  Velocity,  msec 

14.4 

25.6 

11.2 

Average  Acceleration 
from  Initial  Motion 
to  Peak  Velocity,  g 

60 

53.6 

64.8 

In  addition  to  the  gage  data  for  the  south  headwall  above  the  center 
of  the  door,  experimental  data  exist  also  in  the  form  of  motion  picture  film. 
From  these  films  it  was  possible  to  estimate,  although  not  with  great 
precision,  the  velocity  of  the  upper  corners  of  the  doorways  of  the  south 
and  west  igloos.  Although  it  was  also  possible  to  measure  displacements, 
flying  debris  in  the  igloos  obscured  the  view  of  the  door  prior  to  maximum 
displacement.  The  peak  velocities  of  the  corner  points  are  presented  in 
Table  5-2. 


TABLE  5-2.  PEAK  VELOCITIES  OF  UPPER  CORNER  OF  DOOR 


r 

* 

Source 

Peak  Velocity,  fps 

|. 

i 

South  Igloo  (Measured) 

33 

| 

West  Igloo  (Measured) 

19 

[ 

Current  Calculation  (Node  49) 

19-2 

i 

£. 

1 

Prior  Calculation  (Node  31) 

16.6 
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Figures  5-6  and  5“9  show  the  motion  of  points  on  the  interior  of 
the  steel  door.  The  displacements  are  very  large.  However,  the  displacement 
values  shown  in  the  figures  are  probably  typical  of  the  order  of  magnitude 
of  the  actual  door  displacements.  Figure  5-19  (from  Ref.  2)  shows  the  motion 
of  the  point  in  the  prior  calculation  corresponding  to  that  of  Figure  5*6. 
Actual  measurements  from  the  motion  picture  film  of  the  west  igloo  are 
available,  although  visibility  was  obscured  prior  to  maximum  displacement. 

The  maximum  deflections  of  the  center  of  the  door,  as  determined  by  the 
three  independent  means,  are  listed  in  Table  5*3* 


TABLE  5-3.  PEAK  DEFLECTION  OF  CENTER  OF  DOORWAY 


Source 

Displacement,  in. 

West  Igluo  (Measured) 

30* 

Present  Calculation  (Node  73) 

20.6 

Prior  Calculation  (Node  49) 

10.4 

*The  motion  picture  record  terminated  prior  to  maximum 
displacement.  The  value  shown  is  the  largest  observed 
displacement  prior  to  termination  of  the  record. 


Figure  5-12  shows  the  displacement  histories  of  several  points  at 
ground  level  on  the  headwall.  Figure  5-20  (from  Ref.  2)  is  comparable  to 
Figure  5“12a.  From  this  comparison  it  is  evident  that  there  is  a great 
difference  between  the  two  calculations  in  modeling  the  foundation.  The 
foundation  of  the  present  calculation  is  much  more  rigid  than  the  founda- 
tion of  the  prior  calculation.  Accelerometers  located  on  the  floors  of  both 
the  south  and  the  west  igloos  provide  experimental  base  motion  data.  The 
results  are  presented  in  Table  5-4. 
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TABLE  5-4.  PEAK  ACCELERATION  OF  BASE  OF  HEADWALL 


Source 

Peak  Acceleration,  g 

South  Igloo  (Measured) 

6.3 

West  Igloo  (Measured) 

5.5 

Present  Calculation  (Node  104) 

19.8 

Present  Calculation  (Node  108) 

12.9 

Prior  Calculation  (Node  52) 

66.2 

5.2.2  STATIC  HEADWALL  MEASUREMENT 

The  static  headwal 1 measurements  involved  setting  up  survey  monuments 
3 ft  in  front  of  the  igloo  headwails.  The  distances  from  the  monuments  to 
the  headwails  were  measured  at  selected  points  before  the  test.  These  dis- 
tances were  again  measured  after  the  test.  The  net  changes  in  position 
represented  the  permanent  displacement  of  the  headwails.  A gap  of  approxi- 
mately 0.2  ft  was  found  between  the  back  of  the  south  headwal 1 and  the  earth 
cover  at  the  top  left  corner,  and  a similar  space  0.1-ft  wide  was  found  at 
the  top  right  corner.  Th'i  demonstrates  that  the  static  measurements  are  not 
the  maximum  displacements  but  are  the  permanent  movements  cf  the  headwal 1. 

These  permanent  displacements  of  the  headwails  are  shown  in  Figures  5-21 
and  5-22  (Ref.  1).  All  the  measurements  are  subject  to  +0.05  ft  of  error 
because  pretest  measurements  showed  that  the  walls  deviated  from  a true 
vertical  plane  by  that  amount.  The  measured  displacement  patterns  indicate 
that  the  headwails  appear  to  have  responded  in  different  ways.  More  pronounced 
yielding  on  the  steel  arch  was  found  in  the  south  igloo.  In  the  west  igloo, 
there  was  clear  indication  that  the  steel  arch  acted  as  a reaction  line 
resisting  headwal 1 movement.  The  difference  in  the  observed  data,  therefore, 
implies  that  material  properties  in  the  two  igloos  are  not  the  same.  As  was 
mentioned  before,  the  fin’te  element  models  used  to  predict  the  -esponses  of 
the  south  and  west  igloos  are  identicai.  The  assumed  material  properties  of 
the  model  would  therefore  seem  to  represent  characteristics  that  are  a 
compromise  between  the  actual  properties  of  the  south  iglco  and  those  of 
the  west  igloo. 
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Figure  5-23  shows  a contour  plot  of  the  computed  displacements  at 
0.22  sec  of  the  current  ESKIMO  I calculation.  As  can  be  seen  from  the 
response  curves  (Figs,  5- I through  5-12),  motion  of  the  headwall  has 
essentially  stopped  by  that  time'.  This  figure  can  be  interpreted,  therefore, 
as  a plot  of  the  computed  permanent  displacements.  The  deformed  shape  of 
the  hiadwalt  appears  to  be  very  similar  to  that  of  the  west  igloo,  although 
the  .nagnitudes  of  the  displacements  differ  considerably.  It  appears  as  if 
the  west  igloo  experienced  a rigid  body  displacement  not  experienced  in  the 
calculation.  The  backfill  and  the  foundation  of  the  finite  element  model 
were  too  stiff  to  permit  substantial  deformation.  If  a rigid  body  displace- 
ment of  0.15  ft  is  added  to  every  point  of  Figure  5-23,  the  contour  plot 
of  Figure  5*2k  results.  This  plot  is  in  much  better  agreement  with  Figure  5-22 
than  is  Figure  5-23.  The  rigid  body  displacement  does  not  change  the  stress 
field  in  the  headwall.  The  largest  discrepancy  between  Figures  5-2^  and  5-22 
is  the  magnitude  of  the  displacements  above  the  door.  The  calculation  shows 
them  to  be  substantially  higher  than  the  measured  values.  This  may  be  due 

to  the  inability  of  hinges  to  fail  in  the  finite  element  model. 

• 

There  seems  to  be  no  means  of  adjusting  Figure  5-23  such  that  it 
compares  well  with  the  deformation  of  the  south  headwall,  shown  in  Figure  5-21. 
It  is  surprising  that  the  south  and  the  west  igloos,  being  of  essentially  the 
same  construction  and  subject  to  similar  load  conditions,  should  exhibit  such 
dissimilar  deformation.  Figure  5”21  indicates  that  the  most  pronounced 
feature  of  the  deformation  of  the  south  headwall  seems  to  be  the  folding 
about  a vertical  axis.  In  the  case  of  the  west  headwall,  there  seems  to 
have  been  a tendency  toward  folding  about  a horizontal  axis.  The  difference 
between  the  computed  deformations  of  the  finite  element  analysis  and  those 
of  the  west  headwaP  appears  to  be  smaller  than  the  difference  between  the 
deformations  of  the  south  and  west  headwalls.  This  implies  that  the  stiffness 
properties  of  the  igloos,  particularly  of  the  soil  backfill,  are  variable. 
Furthermore,  these  variable  properties  exert  considerable  influence  on  the 
deformation  of  the  headwalls. 
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The  fine)  displacements  of  the  current  calculation  are  considerably 
smaller  than  the  peik  displacements.  It  has  already  been  noted  that  strong 
evidence  exists  to  indicate  that  the  same  Is  true  of  the  actual  head-  nils, 
for  the  sake  of  comparison  with  the  final  deformations,  Figure  5-25  shows 
the  displacement  contours  of  the  headwall  of  the  finite  element  model  at 
0.042  sec,  a time  at  which  much  of  the  headwall  was  in  the  vicinity  of  peak 
lisplacement.  The  comparison  indicates  that  the  final  displacement  contours 
do  not  provide  an  accurate  measure  of  the  severity  of  loading  that  occurred 
during  the  blast.  Time~dependent  measurements  are  essential  for  determining 
peak  stresses  during  such  an  experiment. 

5.2.3  CONCLUSIONS 

The  finite  element  calculation  dis  ■ sed  in  this  section  is  a 
modification  of  a previous  calculation  of  the  response  of  the  south  and 
west  headwatls  of  the  ESKIMO  I event.  The  modifications  included  increasing 
the  number  of  elements,  representing  the  soil  backfill  and  the  foundation  by 
different  models,  changing  the  boundary  condition  along  the  bottom  edge  of 
the  door,  and  increasing  the  flexibility  of  the  concrete.  Several  differences 
between  the  present  calculation  and  the  earlier  calculation  have  been  pointed 
out.  The  most  notable  difference  is  the  tendency,  in  the  present  calculation, 
for  final  displacements  to  be  significantly  smaller  than  peak  displacements. 
The  differences  between  the  calculations  appear  to  result  primarily  from 
the  changes  in  the  type  of  model  (for  example,  using  nonlinear  springs  for 
the  backfill)  rather  than  from  the  increase  in  the  number  of  elements.  It 
appears  that  the  finite  element  grids  are  fine  enough,  n.td  that  further 
improvements  in  correlation  between  calculation  and  experiment  are  to  be 
obtained  through  improving  the  means  of  modeling  the  influence  of  the  various 
components . 
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Comparison  of  the  calculated  results  with  the  experimental  resu'ts 
indicates  that  improved  correlation  might  be  obtainable  by  modeling  breakable 
hinges  along  the  top  and  side  boundaries  of  the  door.  This  would  limit  the 
influence  of  the  door  on  the  headwatl,  and  perhaps  lower  the  headwel!  dis- 
placements adjacent  to  the  doorway,  in  order  to  incorporate  such  a hinge 
element  into  future  calculations,  it  will  be  necessary  to  know  the  ultimate 
strength  of  the  door/headwal 1 connection. 

A second  area  of  modification  of  the  finite  element  model  is  the 
stiffness  of  the  soil  backfill  and  the  foundation.  The  permanent  displace- 
ments of  the  headwails  indicate  that  these  components  should  be  modeled  as 
softer.  As  pointed  out  in  Section  4,  assigning  values  to  the  soil  parameters 
was  largely  a matter  of  intuition,  plus  trial  and  error,  experimental 
measurements  establishing  the  nature  of  participation  of  the  soil  in  the 
response  of  the  headwall  are  needed.  This  participation  is  evidently  quite 
complex,  thus  rendering  static  laboratory  tests  of  soil  properties  inadequate. 
In  fact,  the  only  experimental  means  of  determining  che  participation  of  the 
backfill  seems  to  be  the  ESKIMO  series  itself.  It  must  be  noted  that  further 
reduction  of  the  stiffness  of  the  backfill  will  result  in  the  steel  arch 
applying  a significant  influence  on  the  headwall.  Thus  future  calculations 
may  have  to  Include  the  effects  of  the  arch. 

The  comparison  of  final  displacements  of  the  south  and  west  headwails 
indicated  a variability  in  the  stiffness  of  various  components  of  the  igloos. 
This  variability  increases  the  confusion  over  the  role  of  the  backfill,  the 
arch,  ..«d  the  foundation  in  the  headwall  response.  It  suggests  that  there 
may  always  be  a significant  degree  of  uncertainty  concerning  the  parameters 
used  in  modeling  a headwall  test. 
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5.3  RESPONSE  OF  EAST  HEADUAll,  ESKIMO  II 

The  finite  element  model  for  the  eest  headwall  of  ESKIHO  II  is  the 
seme  es  that  used  for  the  south  heedwell  of  ESKIHO  I , since  the  two  heedwells 
ere  of  identicel  construction.  Therefore,  the  model  shown  in  Figure  3-2 
(Sec.  3)  wes  subjected  to  the  pressure  loading  applicable  to  the  ESKIHO  II 
test.  The  pressure  loading  histories  are  shown  in  Figure  2-3  (Sec.  2).  The 
results  from  the  dynamic,  response  calculation  are  presented  as  time-history 
plots  of  the  headwatl  motion  and  contour  plots  of  the  displacement  of  the 
headwall.  Comparisons  are  made  between  the  calculation  and  the  available 
experimental  measurements  on  the  east  headwall  in  ESKIHO  II. 

5.3.1  RESPONSE  TIHE  HISTORIES 

Since  the  model  used  in  this  study  is  the  same  as  the  one  used  for 
the  south  headwall  of  ESKIHO  I except  for  the  pressure  loading,  the  response 
time  histories  are  shown  at  the  same  nodes  as  in  Section  5.2  in  order  to 
compare  the  results  due  to  different  loading  environments.  These  responses 
are  plotted  in  Figures  5-26  through  5“36. 

The  general  shapes  of  these  time  histories  are  similar  to  the 
corresponding  ones  for  .the  south  headwall  in  ESKIMO  I except  that  the 
magnitudes  of  the  present  time  histories  are,  in  most  cases,  greater.  It 
is  seen  from  Figure  5-26,  for  example,  that  the  displacement  of  Node  1 rises 
initially  to  a positive  peak  of  0.8  in.,  then  reverses  the  motion  to  reach 
a negative  peak  of  2.5  in.  The  positive  peak  is  25%  greater  than  the  corre- 
sponding peak  at  the  same  node  for  the  south  headwall  in  ESKIMO  I,  whereas 
the  negative  peak  is  k0%  greater. 

The  motions  at  nodes  located  outside  the  steel  arch,  where  the 
headwail  is  supported  by  soil,  are  only  slightly  greater  than  those  in 
ESKIMO  I.  However,  the  peak  pressure  In  ESKIMO  II  is  roughly  three  times 
greater  chan  that  in  ESKIHO  I.  Within  the  arch,  the  headwall  motions  in  tne 
present  case  are  much  greater.  For  example,  at  Node  50  th'  rjximum  dis- 
placement is  about  two  and  a half  times  greater  than  the  correspond i ng  value 
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In  ESKIMO  I;  and  at  the  canter  of  the  door  (Node  73).  the  maximum  displacement 
Is  about  three  times  greater.  This  phenomenon  may  be  due  to  the  fact  that  the 
headwall  and  the  door  within  the  arch  are  unsupported  and  yield  sooner  than 
the  region  outside  the  arch.  After  yielding,  additional  loading  causes 
disproportionately  large  displacements. 

Figure  5-33*  shows  that  the  maximum  displacement  of  the  center  of 
the  door  is  about  86  in.,  which  is  greater  than  the  width  of  the  door  (60  in.). 
Since  the  INSLAB  computer  program  is  based  on  the  theory  of  small  displacements, 
the  above  displacement  value  may  be  subjected  to  significant  error.  However, 
the  large  displacements  predicted  by  the  program  simply  reinforce  the  fact  that 
the  door  would  be  bent  out  of  shape  under  the  loading  experienced  by  the  headwall. 

As  noted  above,  the  response  calculation  of  the  east  headwall  showed 
that  the  steel  door  would  be  blown  open  under  the  pressure  loading  assumed 
for  the  ESKIMO  II  test,  After  the  door  is  blown  open,  a redistribution  of 
the  pressure  on  the  headwall  will  take  place,  resulting  in  a net  reduction 
of  the  total  load  on  the  headwall.  However,  the  calculation  assumed  that 
the  pressure  was  acting  uniformly  on  the  headwall  and  the  door  throughout 
the  calculation.  Therefore,  the  calculation  considered  higher  loading  than 
was  the  actual  case  and  the  computed  headwall  notions  would  be  significantly 
different  after  the  door  is  blown  open.  The  computation  for  the  response  of 
the  east  headwall,  ESKIMO  II,  was  terminated  at  about  165  msec;  whereas  for 
the  ESKIMO  I calculation  this  point  was  at  about  2^0  msec.  Because  of  this, 
the  permanent  displacements  of  the  headwall  are  not  clearly  identified  for 
most  nodes  in  Figures  5-26  through  5-36,  In  any  event,  for  the  reason 
mentioned  above,  the  computed  permanent  displacements  are  expected  to  be 
in  significant  error. 
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5.3.2  COMPARISON  OF  COMPUTED  RESULTS  WITH  TEST  DATA 

In  the  ESKIMO  II  test,  linear  motion  transducers  were  positioned 
above  the  center  of  the  doorways  of  the  five  igloos.  For  the  east  igloo, 
this  location  corresponds  to  Node  46  of  the  finite  element  model  (Fig.  3*2). 
A comparison  between  the  test  data  and  the  computed  motions  for  this  loca- 
tion is  shown  in  Table  5*5.  The  computed  duration  and  velocities  are  much 
higher  than  the  measured  values.  However,  as  in  the  case  of  the  south  igloo 
In  ESKIMO  I,  the  average  acceleration  is  in  good  agreement.  As  suggested  in 
Reference  4,  the  measured  data  at  this  location  are  very  sensitive  to  the 
type  of  doorstop  devices  used  for  Inhibiting  the  movement  of  the  door  top, 
relative  to  the  headwall.  Also,  as  discussed  before,  the  computed  motions 
may  be  greater  because  the  separation  of  the  door  from  the  headwall  was  not 
considered  in  the  analysis. 

Two  accelerometers,  one  on  each  door  leaf,  were  mounted  on  the 
east  igloo  (ESKIMO  II).  However,  these  gages  failed  to  perform  properly 
and  therefore  no  useful  data  were  obtained,  as  observed  in  Reference  4. 

The  static  measurements  intended  to  record  the  permanent  headwall 
deformations  were  perforned  as  described  in  Reference  4.  Figure  5-37  shows 
these  deformations  on  one  half  of  the  headwall  and  door  system  for  the  east 
igloo,  ESKIMO  I I . 

Since  the  measurements  are  made  on  the  whole  headwall  but  only 
one  half  of  the  headwall  is  used  in  the  finite  element  model,  there  are  two 
measured  values  for  each  node  of  the  model  except  those  at  the  centerline 
of  the  headwall.  The  contours  in  Figure  5*37  are  based  on  the  average  of 
the  two  measured  values  where  applicable.  All  the  values  are  in  hundredths 
of  a foot,  and  the  positive  value  indicates  deformation  into  the  igloo. 
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TABLE  5-5.  MOTION  OF  EAST  HEADWALL  (ESKIMO  II) 
ABOVE  CENTER  OF  DOORWAY 


1 tern 

Measured 

Value 

Computed 
Value 
(Node  46) 

Maximum  velocity, 
ft/sec 

33.3 

102 

Average  velocity 
• from  initial  motion 

to  peak,  f t/sec  ‘ 

31.5 

72 

Time  from  initiation 
of  motion  to  maximum 
velocity,  msec 

6.75 

20 

Average  acceleration 
from  initial  motion  to 
maximum  veloci ty , g 

153 

155 
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As  was  discussed  In  Section  5*3t  the  east  headwail,  ESKIMO  II,  had 
not  reached  steady  state  at  the  end  of  the  calculation.  Therefore,  It  Is 
difficult  to  compare  quantitatively  the  calculated  results  with  measured 
results.  However,  all  the  nodes  indicate  that  the  accelerations  and  velocities 
become  small  near  the  end  of  the  computation.  The  displacement  contours  at 
two  instants  of  time  ’ate  In  the  calculation  (t  ■ 145  msec  and  t - 157.5  msec) 
are  shown  in  Figures  5*30  and  5“ 39 - In  both  figures  the  maximum  displacements 
are  at  the  top  center  of  the  door,  with  an  order  of  magnitude  of  1.3  ft. 
Although  this  is  much  higher  than  the  measured  value,  the  general  deformed 
patterns  are  quite  similar  to  the  measurements  (Fig.  5*37).  In  addition, 
the  contours  in  Figure  5*37  do  not  extend  to  the  doorway,  because  the  door 
was  blown  off  during  the  test.  The  computed  displacement  contours,  however, 
extend  routinely  into  the  doorway. 

5.4  RESPONSE  OF  NORTHEAST  HEADWALL.  ESKIMO  II 


This  section  presents  the  results  from  the  dynamic  analysis  of 
the  northeast  headwal!  in  ESKIMO  II.  The  results  from  the  dynamic  response 
calculation  are  presented  as  time-history  plots  of  the  headwail  motion  and 
contour  plots  of  the  displacement  of  the  headwail.  Also  presented  are 
comparisons  of  experimental  data  with  the  results  from  the  calculation.  The 
finite  element  model  of  the  headwail  is  shown  in  Figure  3*16  (Sec.  3).  This 
model  is  subjected  to  the  blast  pressure  loadings  shown  in  Figure  2-3 
(Sec.  2).  These  same  loadings  are  applicable  to  the  east  headwail,  since 
the  two  headwal Is  were  at  the  same  distance  from  the  source. 


5.4.1  RESPONSE  TIME  HISTORIES 

Figures  5-40  through  5*55  show  displacement,  velocity,  and  acceiera 
tion  time  histories  computed  at  several  nodes  in  the  finite  element  model 
of  the  northeast  headwail  and  door  system  of  ESKIMO  !!  As  in  the  case  of 
the  south  headwail,  ESKIMO  I,  the  displacement  histories  for  this  case 
exhibit  lower  frequencies  compared  to  the  input  loading.  The  rise  times  of 
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the  displacement  histories  are  longer  than  20  msec,  whereas  the  rise  time  of 
the  input  pressure  history  is  1 msec.  As  observed  in  Section  5.2.1,  this 
phenonenon  is  attributed  to  the  low  fundamental  frequency  of  the  headwall  and 
door  sys  :em. 


I 
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Table  5-6  shows  tne  computed  maximum  and  permanent  displacements  at 
several  nodes.  It  is  clear  from  this  table  that  the  maximum  displacements  of 
the  nodes  inside  the  arch  line  are  much  greater  than  the  corresponding  values 
of  the  nodes  outside  the  archline.  However,  no  such  correlation  exists  for 
the  permanent  displacements.  The  observation  implies  that  the  soil  behind 
the  headwall  Is  very  effective  in  resisting  the  dynamic  forces  on  the  wall, 
but  does  not  greatly  influence  the  residual  displacements  of  the  wall. 

Although  the  east  headwall  is  subjected  to  the  same  loading  as  the 
northeast  headwall,  the  geometries  of  these  two  headwalls  are  different. 
Therefore,  the  relative  responses  of  the  headwalls  will  be  compared  in  order 
to  get  an  irsight  into  their  relative  strengths  and  their  performances  in 
the  test.  Because  of  the  differences  in  the  finite  element  models  of  the 
headwalls,  a node  in  one  model  may  not  correspond  exactly  to  a node  in  the 
other.  Each  comparison  is  therefore  made  between  two  nodes  located  in  the 
same  general  area. 

For  the  northeast  headwall.  Node  9 Is  located  at  the  crown  of  the 
steel  arch.  The  maximum  displacement  at  this  node  is  4.7  in.  The  corre- 
sponding node  for  the  east  headwall,  Node  19,  has  the  maximum  displacement 
of  8.2  in.  Similar  comparison  is  seen  in  Table  5-7  for  other  locations  on 
or  inside  the  arch.  Although  the  area  of  the  northeast  headwall  supported 
by  the  backfill  soil  is  much  smaller  than  that  of  the  east  headwall,  the 
additional  stiffness  due  to  the  presence  of  the  concrete  beams  around  the 
doorway  of  the  northeast  headwall  appears  to  have  reduced  tne  maximum 
displacements  of  the  headwall  significantly. 
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TABLE  5-6.  COMPUTED  MAXIMUM  AND  PERMANENT  DISPLACEMENT 
OF  NORTHEAST  HEADWALL,  ESKIMO  II 


Node  No. 


Max I mum 
Displacement, 


Permanent 
Di splacement, 
in. 


Location 


On  the  arch  line 

Outside  the  arch 

On  the  arch  line 

Inside  tne  arch, 
on  door  frame 

Inside  the  arch, 
on  door  frame 

On  the  arch  line 

Outside  the  arch 

Inside  the  arch, 
on  center  of  door 

Inside  the  arch, 
on  door  frame 

Inside  the  arch 

On  the  arch  1 I ne 

Outside  the  arch 

On  the  arch  I ine 

Inside  the  arch, 
on  door  frame 

Inside  the  arch, 
on  door  frame 

Outside  the  arch 
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TABLE  5-7.  COMPARISON  OF  COMPUTED  MAXIMUM 
DISPLACEMENTS  FOR  ESKIMO  II 


Eas  t Headwa 1 1 

Northeast  Headwail 

n 

Node 

Maximum 

Displacement, 

in. 

Node 

Max i mum 
Di  splacement , 
in. 

General  Location 

n 

8.2 

9 

4.7 

Crown  of  the  arch 

i 

45. C 

39 

7.3 

Top  center  of  door 

49 

20.0 

42 

6.0 

Top  corner  of  door 

73 

86.0 

72 

17.0 

Center  of  door  • 

76 

21.0 

75 

4.7 

Center  edge  of  door 

77 

15.0 

•77 

4.0 

Inside  arch,  along 
horizontal  centerline 
of  door 

103 

0.3 

108 

0.4 

Bottom  corner  of  door 
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Perhaps  the  most  obvious  difference  between  the  two  headwali 
responses  is  seen  in  the  steel  door  area.  The  maximum  displacements  in 
the  door  area  for  the  northeast  headwali  are  much  smaller  than  those  for 
the  east  headwali  (Table  5*7).  At  the  center  of  the  door,  for  example,  the 
maximum  displacement  for  the  northeast  igloo  is  17  in.  as  against  86  in. 
for  the  east  headwali.  As  described  in  Section  1,  the  door  system  for  the 
northeast  igloo  is  of  the  biparting  and  sliding  type,  whereas  the  door 
system  for  the  east  igloo  is  of  the  double  leaf  and  hinged  type.  The  results 
of  the  calculations,  therefore,  suggest  that  the  biparting  and  sliding  type 
is  superior  to  the  double  leaf  and  hinged  type  for  resisting  the  blast  loads. 

5. A. 2 COMPARISON  OF  COMPUTED  RESULTS  WITH  TEST  DATA 

As  explained  in  Section  5*3*2,  linear  motion  transducers  were 
placed  above  the  center  of  the  doorways  of  all  the  igloos  in  the  ESKIMO  U 
test.  For  the  northeast  igloo  this  position  corresponds  to  Node  39  of  the 
finite  element  model  (Fig.  3“ 1 6 , Sec.  3).  A comparison  between  the  test 
data  and  the  computed  motions  for  this  location  is  shown  in  Table  5-8. 

Among  the  four  items  compared,  the  times  from  initiation  of  motion  to  maximum 
velocity  are  in  good  agreement  between  the  computed  and  measured  values.  The 
remaining  items  show  that  the  computed  values  are  abouc  two  times  greater. 

As  mentioned  before,  the  motion  of  this  point  is  very  sensitive  to  the  type 
of  door-stop  device  used  for  inhibiting  the  movement  of  the  door  relative 
to  the  headwali.  Also,  the  computed  motions  may  be  greater  because  the 
separation  of  the  door  from  the  headwali  was  not  considered  in  the  analysis. 

Two  accelerometers  were  mounted,  one  on  each  door  leaf.  The  gage 
on  the  left  door  leaf  failed  to  perform,  while  the  gage  on  the  right  door 
leaf  recorded  a peak  acceleration  value  of  1200g.  The  computed  value  at 
the  corresponding  location  (Node  72)  is  seen  from  Figure  5-^7  to  be  1300g. 
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TABLE  5-8.  MOTION  OF  NORTHEAST  HEADWALL,  ESKIMO  II, 
ABOVE  CENTER  OF  DOORWAY 


Max i mum  ve  < oc i ty , 
ft/sec 

Average  velocity 
from  initial  motion 
tc  peak,  ft/sec 

Time  from  initiation 
of  motion  to  maximum 
velocity,  msec 

Average  acceleration 
from  initial  motion  to 
maximum  velocity,  g 


Measured 

Value 

Computed 
Value 
(Node  39) 

14.3 

32 

10.4 

21 

17 

19 

26 

52 

The  permanent  deformations  of  the  northeast  headwall  are  given  in 
Reference  4,  and  are  shown  plotted  in  Figure  5~56.  The  contours  in  the 
figure  are  based  on  the  average  of  the  twc  measurements  at  symmetrical 
points  on  the  headwall.  Figure  5 “57  shows  the  computed  displacement 
contours  of  the  headwall  at  109  msec.  Since  the  motions  become  essentially 
steady  state  at  this  time,  the  displacements  of  the  headwall  at  this  instant 
are  equal  to  the  permanent  displacements  of  the  headwall. 


Figures  5“56  and  5-57  show  that  the  maximum  deformations  occur 
around  the  doorway.  The  contours  seem  to  follow  the  patterns  of  the  steel 
arch  and  the  floor  slab.  In  general,  the  measured  values  are  slightly 
higher  than  the  computed  values.  In  both  figures,  the  movements  tend  to  be 
smaller  toward  the  top  and  the  bottom  of  the  headwall.  At  the  top  of  the 
headwall,  the  computed  displacements  are  negative,  whereas  the  measurements 
show  small  positive  values.  In  general,  the  wingwall  experiences  smaller 
displacements.  At  the  far  end  of  the  wingwall,  some  displacements  are 
negative.  This  may  be  due  to  the  fact  that  the  headwall  is  supported  by 
the  soil  outside  the  arch  and  is  subject  to  rotation  about  the  support. 


5-5  RESPONSE  OF  NORTHEAST  HEADWALL,  ESKIMO  IV 

The  response  calculation  of  the  northeast  headwall,  ESKIMO  IV,  was 
to  have  been  performed  before  tne  test  was  conducted  on  10  September  1975, 
so  that  the  analytical  results  could  be  used  to  predict  the  behavior  of  the 
headwall  and  door  system  in  the  test.  Because  of  tight  schedules,  the  above 
calculation  could  not  be  completed  before  the  test;  but  the  calculation  was 
performed  before  any  results  from  the  test  were  available.  This  assured 
that  the  calculation  was  performed  independently  of  the  test  and  was  not 
subject  to  any  bias  that  could  have  resulted  from  reviewing  the  test  results 
before  performing  the  calculation. 
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As  of  the  writing  of  this  report,  the  data  from  the  ESKIMO  IV  test 
ere  evei table  only  In  preliminary  form.  Therefore,  the  test  results  reflect 
raw  data  that  may  be  subject  to  future  corrections.  This  fact  should  be 
kept  In  mind  when  comparing  the  analytical  results  with  the  test  data. 

5.5.1  RESPONSE  TIME  HISTORIES 

Figures  5—58  through  5*69  show  displacement,  velocity,  nd  accelera- 
tion time  histories  computed  at  several  locations  on  the  northeast  headwall 
of  the  ESKIMO  IV  event.  These  locations  generally  correspond  to  those  of  the 
northeast  headwall  of  the  ESKIMO  II  event,  for  which  the  motion  time  histories 
are  shown  in  Figures  5*40  through  5*55.  Since  the  structural  configuration 
of  the  no*theast  headwall,  ESKIMO  IV,  closely  resembles  the  northeast  headwall 
of  the  ESKIMO  II  test,  a comparison  of  the  responses  of  both  northeast 
headwa Us  will  be  made . 

The  motion  time  histories  shown  in  Figures  5*58  through  5-69  reflect 
calculations  with  a duration  of  1 86  msec.  The  displacement  time  histories  at 
most  locations  show  a single  peak  at  about  25  msec  and  become  steady  state 
after  about  100  msec.  These  time  histories  are  quite  similar  to  those  for 
the  northeast  headwall  of  ESKIMO  II,  both  in  magnitudes  and  general  behavior. 
The  input  blast  pressure  for  ESKIMO  IV  is  slightly  higher  than  that  for 
ESKIMO  II,  but  the  resulting  peaK'di splacements  of  the  headwall  sre  slightly 
lower  in  ESKIMO  IV,  as  shown  in  Table  5*9.  Therefore,  it  can  be  concluded 
that  the  northeast  headwall  structure  in  ESKIMO  IV  is  stronger  than  that  in 
ESKIMO  II  for  resisting  blast  loadings. 

The  response  time  histories  correspond! ng  to  the  steel  door  at 
Nodes  61,  64,  91,  94,  121,  and  124  are  shown  respectively  in  Figures  5*60, 

5-61 , 5—63 , 5*64,  5-68,  and  5*69.  For  the  northeast  igloo  of  ESKIMO  II, 
the  oeak  displacements  along  the  edges  of  the  door  are  only  slightly  higher; 
but  those  at  the  center  of  the  door  are  significantly  highe’r  than  the  values 
shown  in  Table  5*9  for  ESKIMO  IV.  The  primary  reason  for  these  differences 
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TABLE  5-9-  COMPARISON  OF  MAXIMUM  DISPLACEMENTS  OF  NORTHEAST 
HEADUALLS  OF  ESKIMO  II  AND  ESKIMO  IV  IGLOOS 


ESKIMO  II 

ESKIMO  IV 

Node 

Max i mum 
Displacement, 
in. 

Node 

1 

Max i mum 
Displacement, 
in. 
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is  the  door  configurations.  In  ESKIMO  IV,  th«  door  i s of  a single-leaf 
sliding  typa,  whrraas  tha  door  in  ESKIMO  II  is  of  a biparting  sliding  typa. 

Tha  single-leaf  typa  has  a much  graatar  sciffnass  tc  rasist  daf taction  at  tha 
cancar  of  tha  door  than  does  tha  doubte-leaf  typa.  Tharafora,  othar  things 
baing  aqua),  tha  double- leaf  door  experiences  highar  daf  I action  than  tha 
single-laaf  ona. 

5.5.2  COMPARISON  OF  COMPUTED  RESULTS  WITH  TEST  DATA 

As  in  tha  casa  of  tha  othar  Eskimo  tasts,  posttast  measurements 
were  mada  on  tha  haadwal Is  of  tha  ESKIMO  IV  tast  (Raf.  6),  in  ordar  to 
obtain  parmanant  dispiacamants  of  tha  haadwal Is.  Figura  5*70  shows  tha 
parmanant  displacamant  contour  linas  of  tha  northaast  haadwal I basad  on  tha 
avaraga  of  two  maasurad  valuas  dua  to  tha  symmatry  of  tha  haadwal I at  tha 
vartical  cantarlina.  Tha  camputad  displacamant  contours  dt  tima  t * 158  msec 
ara  shown  in  Figura  5-7l . As  notad  bafora,  tha  computad  dispiacamants  become 
staady  stata  aftar  about  100  msac  and  tharafora  tha  dispiacamants  of  tha 
haadwal 1 at  158  msac  raprasant  tha  parmanant  dispiacamants.  A comparison 
of  tha  computed  and  maasurad  contours  shows  that  the  computed  displacements 
generally  ara  smaller  than  those  measured.  However,  the  differences  are  not 
great  and  are  not  significant  because  the  magnitudes  of  the  permanent  displace- 
ments ara  small.  Tha  contours  do  not  indicate  any  strong  influence  of  the 
arch  on  the  haadwall  response. 

Figures  4-3  and  4-4  (Sec.  4)  show  the  locations  of  the  accelero- 
meters and  the  linear  motion  transducers  mounted  on  the  northeast  headwal! 
and  door  for  ESKIMO  IV.  The  measured  acceleration  and  displacement  tire 
histories  for  tha  northaast  haadwall  are  shown  respectively  in  Figures  5-72 
and  5*73.  These  records  ara  reproduced  from  the  preliminary  report  (Ref.  6) 
on  the  ESKIMO  IV  tast.  As  such,  these  records  have  not  been  corrected  for 
drifts  and  othar  random  errors.  This  fact  should  be  kept  in  mind  in  the 

v e 

following  discussion  involving  the  measured  time  histories. 
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The  measured  displacement  time  histories  at  several  locations  on 
the  headwall  anc  door  are  shown  in  Figure  5 -73 - The  displacements  shown  are 
in  inch  units  and  nre  plotted  as  functions  of  time  expressed  in  msec  units. 

Next  to  each  time  history  the  node  in  the  finite  element  mesh  of  the  headwall 
(Fig.  3- 1 8 , Sec.  3)  to  which  the  time  history  corresponds,  is  shown.  Of  the 
eight  measured  displacement  time  histories,  direct  comparison  with  the 
calculated  histories  is  possible  at  only  four  locations.  These  are  at  Nodes  21, 
64,  94,  and  97  and  the  corresponding  time  histories  are  shown  respectively  in 
Figures  5-58,  5-61,  5-64,  and  5-65.  A comparison  of  the  calculated  and 
measured  displacement  at  the  four  locations  shows  that  shapes  of  the  pulses 
are  in  good  agreement.  Table  5-10  shows  a comparison  of  the  maximum  displace- 
ments and  permanent  displacements.  The  agreement  between  the  computed  and 
measured  permanent  displacements  is  very  good.  The  computed  maximum  dis- 
placements are  somewhat  higher  than  those  measured.  Nevertheless,  the 
differences  between  the  computed  and  measured  values  are  small  in  comparison 
with  those  found  for  the  headwalls  in  the  ESKIMO  I and  II  tests. 

The  measured  acceleration  time  histories  at  several  locations  on 
the  headwall  and  door  are  shown  in  Figure  5-72.  The  accelerations  are  given 
in  the  units  of  g (gravity  constr.it)  and  the  time  is  expressed  in  msec  units. 
Next  to  each  record  the  corresponding  node  in  the  Finite  element  mesh 

(Fig.  3- 1 8 , Sec.  3)  is  shown.  Out  of  a total  of  eleven  records,  the  computed 

pulses  are  available  only  at  seven  locations.  These  are  at  Nodes  21,  6 1 , 64, 

91,  94,  97,  and  99  and  the  corresponding  acceleration  pulses  are  given 
respectively  in  Figures  5-58.  5-60,  5-61 , and  5-63  through  5 “66 . The 
measured  pulses  are  shown  up  to  a maximum  of  50  msec,  while  the  duration 
of  the  computed  time  histories  is  1 86  msec.  The  general  shapes  of  the 
computed  pulses  are  different  from  those  measured.  The  measured  records 
show  small  motions  after  about  25  msec.  In  contrast,  the  computed  accelera- 
tion records  show  strong  motions  even  after  150  msec.  This  is  attributed 
to  the  fact  that  in  the  actual  test,  cracking  and  crushing  of  concrete 

dissipates  large  amounts  of  energy.  Although  the  INSLAB  code  accounts  for 
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TABLE  5-10.  COMPARISON  OF  MAXIMUM  AND  PERMANENT  DISPLACEMENTS 
OF  NORTHEAST  IGLOO,  ESKIMO  IV 


Location 

Max ! mum 
Displacements , 
In. 

Permanent 
Displacements , 
in. 

Node 

Gage  1.0. 

Measured 

Computed 

Measured 

Computed 

94 

L-1-B 

3.0 

3.7 

0.8 

mm 

97 

L-2-8 

3.2 

3.3 

1 .0 

HB 

64 

L-5-B 

3.6 

5.3 

1.0 

21 

L-9-B 

1.9 

3.1 

* 

BB 

*The  transducer  seemed  to  have  failed  after  recording  about  50  msec 
of  motion. 
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the  energy  dissipation  tnrcugh  the  yielding  of  concrete.  It  is  not  intended 
to  consider  the  crushing  of  concrete.  Therefore,  the  energy  dissipation 
mode)  in  the  calculation  may  underestimate  the  actual  energy  loss.  A 
comparison  of  the  computed  and  measured  maximum  accelerations  at  various 
locations  on  the  headwall  and  door  Is  shown  in  Table  5-11.  The  agreement 
between  the  measured  and  computed  values  is  good  and  the  differences  *.*e 
small  in  comparison  with  similar  differences  in  the  ESKIMO  I and  II  tes::. 

In  ESKIMO  IV.  unlika  In  the  previous  ESKIMO  tests,  the  headwalls 
were  heavily  instrumented  to  measure  pressure,  acceleration,  and  displacement 
time  histories  on  the  headwall  and  door  systems.  As  described  in  Section  A, 
the  northeast  headwall  was  instrumented  with  seven  pressure  gages,  twelve 
accelerometers , end  nine  linear  motion  transducers.  This  provided  for  the 
first  time  an  c.pportuni  ty  to  compare  the  computed  and  measured  time  histories 
of  accelerations  and  displacements  at  several  locations  on  the  headwall. 

The  information  contained  in  the  measured  records  can  be  used  to  adjust  the 
parameters  of  sci 1 models  in  any  future  calculations  to  improve  the 
analytical  prediction  capabilities. 


A 


R-7556-1-4182 


FIGURE  5-23.  COMPUTED  DISPLACEMENTS  (IN  HUNDREDTHS  OF  A FOOT)  OF  THE  HEADWAILS 
OF  THE  SOUTH  ANO  WEST  IGLOOS,  ESKIMO  I,  0.22  SEC  AFTER  ARRIVAL 
OF  THE  BLAST  WAVE 
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FIGURE  5-24.  COMPUTED  DISPLACEMENTS  (IN  HUNDREDTHS  OF  A FOOT)  WITH  A RIGID 
BODY  DISPLACEMENT  OF  0.15  FT  ADDED  IN  FOR  THE  SAKE  OF 
COMPARISON  WITH  THE  WEST  IGLOO,  ESKIMO  I 
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(a)  Displacement 


(b)  Velocity 


(c)  Acceleration 


FIGURE  5-36. 


COMPUTED  MOTION  OF  EAST  HEADWALL  (ESKIMO  II)  AT  NODE  103 
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FIGURE  5-38.  COMPUTED  DISPLACEMENT  CONTOURS  OF  EAST  HEADWALL  AT  t « 145  MSEC 
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FIGURE  5-39.  COMPUTED  DISPLACEMENT  CONTOURS  OF  EAST  HEADWALL  AT  t - 157  MSEC 
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FIGURE  5-57.  COMPUTED  PERMANENT  DISPLACEMENTS  OF  NORTHEAST  HEADWAI.L  OF  ESKIMO  III 
(IN  HUNDREDTHS  OF  FEET) 
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FIGURE  5-70.  MEASURED  PERMANENT  DISPLACEMENT  CONTOURS  OF  NORTHEAST 
HEADWALL  OF  ESKIMO  IV  (IN  HUNDREDTHS  OF  FEET) 
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CORRESPONDS  TO  NODE  94 
(FIGURE  5-64) 


CORRESPONDS  TO  NOOE  97 
(FIGURE  5-65) 
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SECTION  6 

CONCLUSIONS  AND  RECOMMENDATIONS 


6.1  CONCLUSIONS 


Th«  present  study  consisted  of  dynamic  analyses  of  magazine  head- 
walls  subjected  to  blast  loadings  in  the  ESKIMO  tests.  Ail  the  analyses 
were  performed  using  the  INSLA6  code,  a nonlinear  finite  element  computer 
program.  The  computed  results  were  compared  with  available  test  data.  The 
following  is  a summary  of  conclusions  of  the  study. 

a.  The  modified  INSLA6  computer  program  is  capable  of  analytically 
predicting  the  behavior  of  the  magazine  headwalls  used  in  the 
ESKIMO  tests. 

b.  The  proposed  nonlinear  material  model  to  represent  the  soil 
behind  the  headwalls  appears  to  be  satisfactory. 

c.  The  general  behavior  of  the  heajlwalls  predicted  by  the  calcula- 
tions agreed  well  with  the  actual  behavior  of  the  headwalls 

. in  the  tests. 

d.  The  results  from  the  present  analysis  of  the  south  and  west 
headwalls  in  ESKIMO  I,  using  the  refined  finite  element  mesh, 
show  better  agreement  w«th  the  test  results  than  those  from 
the  previous  analysis. 

e.  The  analysis  of  the  east  headwall  in  the  ESKIMO  II  test  shows 
that  (1)  the  computed  maximum  motions  are  much  higher  than 
those  measured,  and  (2)  the  computed  permanent  displacement 
contours  are  similar  in  shape  to  those  measured. 


Preceding  page  blank 
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f.  The  results  of  the  analysts  of  the  northeast  headwatl  tn 

the  ESKIMO  II  test  show  that: 

1.  The  computed  maximum  velocity  above  the  doorway  is  much 
higher  than  that  measured  using  the  linear  motion 
transducers. 

2.  The  computed  maximum  acceleration  of  the  door  is  in  good 
agreement  with  that  measured  using  the  accelerometers. 

3.  The  presence  of  the  concrete  beams  around  the  doorway  is 
effective  in  reducing  the  maximum  displacements  of  the 
headwa 1 1 . 

4.  The  biparting  and  sliding  type  of  door  is  superior  to  the 
double- leaf  and  hinged  type,  for  resisting  the  blast 
loads. 


g.  The  results  of  the  analysis  of  the  northeast  headwa 1 1 in  the 

ESKIMO  IV  test  indicated  that: 

1.  The  corrugated  steel  arch  doesn't  seem  to- provide  any 
significant  resistance  to  the  motion  of  the  headwall. 

2.  The  magnitudes  and  general  shapes  of  the  computed  dis- 
placement time  histories  are  in  good  agreement  witn 
those  measured  by  the  linear  motion  transducers. 

3.  The  computed  acceleration  time  histories  look  somewhat 
different  from  those  measured  using  the  accelerometers; 
however,  the  computed  maximum  accelerations  are  in  good 
agreement  with  those  measured. 
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6.2  RECOMMENDATIONS 

The  analytical  prediction  capabilities  exhibited  by  the  present 
study  have  been  satisfactory.  However,  there  are  certain  aspects  to  the 
modeling  of  the  headwall  and  door  system  that  can  improve  the  analytical 
prediction  capabilities.  The  following  is  a summary  of  recommendations  for 
improving  the  analytical  prediction  capabilities  in  order  to  reduce  the 
dependence  on  the  expensive,  full-scale  testing  program. 

a.  The  present  techniques  for  analyzing  the  headwall  systems  are 
not  satisfactory  when  the  doors  are  blown  open.  Therefore, 
the  analytical  methods  should  be  improved  to  include  the 
concept  of  breakable  hinges.  This  would  permit  the  calcula- 
tion of  true  response  of  the  headwall  even  after  the  doors 
are  blow  open. 

b.  The  parameters  in  the  nonlinear  soil  models  should  be  adjusted 
based  on  the  test  date  available  for  the  ESKIMO  IV  test.  The 
resulting  models  are  expected  to  better  represent  the  backfill 
soils  for  improved  analytical  predictions. 

c.  Detailed  analysis  of  the  test  data  for  ESKIMO  IV  and  comparison 
of  the  data  with  the  results  of  the  present  analysis  should 

be  undertaken.  This  may  include  corrections  to  the  data  for 
drifts,  etc.,  and  also  statistical  analyses  of  the  test  data. 

d.  Future  tests  should  include  a sufficient  amount  of  instrumenta- 
tion to  produce  data  necessary  to  study  the  stress  distribution 
in  the  headwall  sections.  This  would  include  placing  strain 
?ages  on  reinforcing  bars. 
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APPENDIX  A 

INCLUSION  OF  NONLINEAR  SPRINGS  IN  INSLAB  CODE 


The  response  of  ths  magazine  headwall  to  blast  loading  is  greatly 
influenced  by  the  earth  fill  behind  the  headwall.  The  effect  of  the  soil  can 
best  be  Incorporated  In  the  two-dimensional  headwall  model  by  using  a com- 
bination of  damping  elements  and  nonlinear  (one-way)  springs.  Linear  springs, 
which  are  already  present  In  the  existing  INSLAB  code,  are  inadequate  for 
modeling  the  influence  of  the  soil  because  they  exert  both  pressures  and 
tensile  forces  on  the  headwall  and  they  tend  toward  a final  displacement  of 
zero.  The  properties  of  the  nonlinear  springs  and  the  details  • their  incor' 
poration  in  INSLAB  are  discussed  herein. 

Figure  A-l  shows  the  load/def lection  characteristics  of  the  non- 
linear springs.  Two  different  spring  constants,  k^  and  k^,  appropriate 
in  the  virgin  loading  zone  (Zone  2)  and  the  unloading/reloading  zone  (Zone  1) 
respectively,  have  been  incorporated  in  such  a way  that  after  loading,  equi- 
librium corresponds  to  a nonvanishing  displacement,  6q.  This  equilibrium 

displacement  is  a function  of  5 , the  largest  compressive  displacement  of 

max 

the  prior  load  history.  Tensile  sprjng  forces  are  eliminated  by  using  a 
spring  constant  of  zero  in  "Zone  0,"  in  which  the  displacement  is  on  the 
tensile  side  of  the  equilibrium  displacement. 

Tht  varying  status  of  the  nonlinear  spring  between  Zones  0,  1, 
and  2 requires  reformulation  of  the  global  stiffness  matrix  during  transition 
f from  one  zone  to  another.  Furthermore,  as  shown  in  Table  A-l,  transition 

?'  between  zones  usually  requires  an  adjustment  to  the  load  vector.  This  is 

i because  the  stiffness  matrix  *s  capable  of  accomodatin'1  only  that  portion  of 

; the  Incremental  spring  force  that  is  proportional  to  the  incremental  displace- 

J.  ment,  A<5.  Note  that  for  a linear  spring  in  which  there  is  never  <jny 

f transition  between  zones,  there  is  no  need  for  a load  vector  adjustment. 
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k - AN  ARBITRARY  STIFFNESS  VALUE  IN  UNLOADING 
u (CONTROLS  THE  PERMANENT  DISPLACEMENT  OF 
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Successful  incorporation  of  tha  spring  force  (Table  A-l)  in  the 
finite  element  model  is  dependent  upon  the  ability  to  predict  upcoming 
changes  in  zone.  This  is  particularly  important  in  the  case  of  rebonding 
(transition  from  Zone  0 to  Zone  1).  In  this  case  the  spring  constant  changes 
from  0 to  kg.  If  this  change  in  stiffness  is  i.ot  anticipated,  calculation 
of  the  incremental  displacement  will  proceed  as  if  the  spring  constant 
remains  0 throughout  the  time  interval.  This  will  result  in  an  excessively 
large  displacement  Into  Zones  1 or  2.  By  the  time  the  spring  force  correspond- 
ing to  this  large  displacement  is  included  in  the  subsequent  step,  an  incre- 
ment of  potential  energy  (spring  deformation)  has  been  erroneously  added  to 
the  dynamic  system.  Cyclic  repetition  of  this  error  causes  unstable  growth 
of  oscillatory  motion. 

Using  a truncated  Taylor  series  expanded  about  time  t,  the  dis- 
placement at  the  end  of  the  time  interval  At  can  he  estimated  according  to 
the  formula 


X(t  + At)  - Xt  + X'  At  + j Xt(At)2  (A-l) 

The  accuracy  of  this  estimate*  should  improve  with  decreasing  At.  This  fore- 
casted future  displacement  is  used  to  predict  whether  any  changes  in  zone  are 
to  be  anticipated.  If  so,  a stiffness  reformulation  occurs  during  the  next 
step,  and  the  appropriate  adjustment  is  incorporated  in  the  load  vector. 

Figure  A-2  shows  the  test  problem  that  was  used  to  check  out  the 
performance  of  the  nonlinear  spring.  Some  sample  results  of  this  test  problem 
are  shown  in  Figure  A-3.  Figure  A-3a  shows  the  time  dependence  of  the  load 
function,  P,  used  in  Case  1.  The  displacement  of  the  end  of  the  beam  for 
this  case  is  shown  in  Figure  A-?b.  The  external  load  (*  2P)  deflects  the 
beam  until  the  nonlinear  springs  supply  an  equilibrating  force.  When  the 

e 

external  load  is  relaxed,  the  nonlinear  spring  relaxes  also.  However,  due 
to  the  high  unloading  spring  constant,  this  relaxation  of  the  spring  occurs 
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with  only  a small  recovery  of  displacement.  When  the  external  load  has 
relaxed  to  the  extent  that  the  elasticity  of  the  beam  overcomes  it,  the  beam 
debonds  from  the  nonlinear  spring  and  undergoes  free  vibration. 

Figure  A-4a  shows  the  time  dependence  of  the  load  function  used  in 
Case  2,  which  was  designed  to  test  rebonding.  The  results  are  shown  in 
Figure  A-4b.  Debonding  occurs  at  about  30  msec  in  this  case,  after  which 
free  vibration  begins.  The  free  vibration  continues  until  the  reapplication 
of  external  load  at  45  msec.  This  externa!  load  and  the  inertia  of  the 
vibrating  beam  cause  the  beam  to  slam  violently  into  the  relaxed  spring  at 
the  displacement  6 (see  Fig.  A-1).  Temporary  equilibrium  is  found  in  the 
vicinity  of  this  displacement  until  the  external  load  grows  large  enough  to 
force  the  beam  into  Zone  2 of  Figure  A-1.  At  this  point  the  stiffness 
reduces  to  the  virgin  loading  value,  allowing  further  displacement  until  a 
new  equilibrium  is  reached. 

The  results  of  these  two  test  cases  indicated  that  the  nonlinear 
spring  is  performing  as  intended.  Both  debonding  and  rebonding  appear  to  be 
working  correctly.  Furthermore,  the  nonzero  equilibrium  displacement  follow- 
ing a load  cycle  has  been  checked  and  found  satisfactory.  It  is  anticipated 
that  these* proper ti es,  when  used  in  conjunction  with  dampers  and/or  concen- 
trated masses,  will  prove  satisfactory  in  modeling  the  effects  of  the  soil 
beh i nd  the  headwa 1 1 . 
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APPENDIX  B 


PROCEDURE  FOR  CALCULATING  YIELD  MOMENTS 
FOR  A REINFORCED-CONCRETE  SECTION 


The  yield  moment  of  a reinforced-concrete  cross  section  is  defined 
as  the  moment  that  makes  the  steel  in  tension  or  the  concrete  in  compression 
reach  its  yieid  limit.  The  reinforcement  provided  in  the  headwall  section 
is  governed  by  the  former  case;  therefore,  only  the  derivation  for  the  under- 
reinforcement case  (steel-in-tension  yield)  is  described  below. 


A typical  cross  section  of  a reinforced-concrete  section  is  shown 

in  Figure  B— 1 , and  the  strain  diagram  is  shown  in  Figure  B-2,  where  A$  and 

are  the  areas  of  the  steel  reinforcement  in  tension  and  compression, 

respectively.  The  strain  level  in  tension  steel  is  e , in  concrete  is  -z  , 

t c 

and  in  compression  steel  is  e'.  The  strains  are  related  by  the  following 
formula: 


d - Kd 

i?a 


1 - K 
K 


d1  - Kd 
Kd 


d * 


(B-1) 


The  stresses  in  steel  are  as  follows: 


1 - K 

a ■ - no  — ■ stress  in  tension  steel 

5 C In 


f -K 

o'  ■ - no  — -r ■ stress  in  compression  steel 

S C R 


(B-2) 


where  n ■ the  ratio  of  modulus  of  elasticity  of  steel  to  concrete. 
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Denoting  the  area  of  steel  as  a fraction  of  the  area  of  the  concrete 
section , i .e.  , 


p " bd 


P' 


V 

bd 


(8-3) 


the  equation  of  equilibrium  becomes 


o.  A + j'  A'  ♦ a Kbd  - 0 

5 3 S 5 4 C 


(B-A) 


This  can  be  rewritten  after  some  arrangement  into  the  form 


X2  + 2(np  + np')  K - (2np  + 2np'  •j-)  ■ 0 


Therefore  the  solution  of  K becomes 


K • ^/(  2np  + 2np*  — — ) + (np  + np 1 ) 2 - (np  + np') 


(B-5) 


The  yield  moment  for  the  underreinforcement  case  is 


A.  (d  - 5^)  ♦ c's  a:  (d1  - ^) 


s s 


or 


M ■ A jdc 
SJ  s 


(B-6) 
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where 


If  the  solution  of  Equation  B-5  is  such  that 
stress  in  compression  steel  is 


then 


T 

s 


2 


o 

s 


The  equation  of  equilibrium  for  this  case  becomes 


a A 
s s 


A ' + -S-  o Kbd  - E e'  A' 


After  some  rearrangement,  and  the  use  of  Equation  B-3,  K becomes 

1 » 

. K ■ ^jl[np  + (2n  - 1 ) p ' j-]  + [np  + (2n  - 1 ) p ' ] 2 


- [np  + (2n  - 1 ) p ' ] 

The  form  of  Equation  B-6  remains  the  same  except  that  for  this  case, 
defined  as 


K 

3 


(2n  - ,)  K-  2.i 


n(l  - KT 


3' 


(B-7) 

the 

(B-8) 

(B-9) 


(B-10) 
j is 


j 


1 


(B-11) 
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APPENDIX  C 

PROCEDURE  FOR  COMPUTING  EQUIVALENT  MOOULI 
FOR  A REINFORCEO  CONCRETE  SECTION 


REINFORCED  CONCRETE  STRUCTURE 


AA8240 


Assume  that  the  displacements  In  the  reinforced  concrete  structure  shown 

♦ 

above  are  in  the  following  form: 


i(x,y,z)  - - z * -W(x)z 


v(x,y,z)  ■ 0 


w(x,y ,z)  • W(x) 


and  the  strains: 


- - — z ■ - W"(x)z 


e ■ e 
zz  yy 


e ■ e 
xz  xy 
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Assuming  that  the  vertical  stress  c » 0 

22 


(A+2G)e  ♦ A(e  ♦ e ) ■ 0 

22  xx  yy 


22 


" A+2G  (cxx  + *yy) 


IT2G  W"{x)2 


and 


v • (w2g>£„  - »*l2 


(A+2G)  - -±r-  t 

A+2G  xx 


(l-v*) 


I exx 


where 


“ — W*1  (x)  2 


_ El 

iV 


Now  the  moment  is  given  by 


h/2 

-7  2d2 

h/2  xx 


bDV",  x rh/2  2 . , v 

-(x)/  z dz  - OW"  (x) 

•'-h/2 
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T 


-U"(x) 


STRAIN  DISTRIBUTION 


Prom  the  above  strain  distribution, 


“"<«>  ' TTETd 


et  « (l-k)dW'(x) 


Moment  of  the  section  Is  then  given  by 


M » Asjdcs  * AsJ’dEset 


M * A jE  (1-k)d2W"(x)  « 0W"{x) 


0 - AJE  (1-k)d2  - — E-V 

5 5 (1-vZ) 


Therefore  the  effective  E for  the  concrete  section  is  given  by 


A JE  (1-k)d2(l-v2) 


